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Abstract: We present a framework for performing a comprehensive analysis of a large 
class of supersymmetric models, including spectrum calculation, dark matter studies and 
collider phenomenology. To this end, the respective model is defined in an easy and straight- 
forward way using the Mathematica package SARAH. SARAH then generates model files 
for CalcHep which can be used with micrOMEGAs as well as model files for WHIZARD 
and 0' Mega. In addition, Fortran source code for SPheno is created which facilitates the 
determination of the particle spectrum using two-loop renormalization group equations 
and one-loop corrections to the masses. As an additional feature, the generated SPheno 
code can write out input files suitable for use with HiggsBounds to apply bounds coming 
from the Higgs searches to the model. Combining all program provides a closed chain from 
model building to phenomenology. 
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1. Introduction 

With the first collisions at the Large Hadron Collider (LHC), a new era of high energy 
physics has started. The LHC is designed to get new and deeper insights into the funda- 
mental principles governing our physical world. It is not only supposed to find the last 
missing particle of the Standard Model (SM) of particle physics, the Higgs boson, but also 
to discover possible physics beyond the SM. 

Despite having proven itself as a successful and precise description of all experiments 
in particle physics for last 30 years, the SM suffers from several theoretical shortcomings, 
and therefore a host of extensions of the model which lead to such new effects have been 
devised over the last years. Among the most prominent examples of such conceptual flaws 
which plague the model are the hierarchy problem [1, 2] and the lack of a candidate for 
Dark Matter. 

Supersymmetry (SUSY) is arguably among the best studied extensions of the SM [3, 4, 
5] which have been developed to cure these problems. Already the minimal, supersymmetric 
extension of the SM, the MSSM (Minimal Supersymmetric Standard Model) [6, 7, 8], solves 
the hierarchy problem, is capable of providing a dark matter candidate [9, 10], leads to 
gauge coupling unification [11, 12, 13, 14] and relates electroweak symmetry breaking to 
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the large top mass [15, 16]. Because of these appealing features, the MSSM was extensively 
studied in the last decades and many software tools were developed for that purpose. 

However, the MSSM might also not be the final answer. One reason is that the MSSM 
like the SM can't explain neutrino masses [17, 18, 19, 20]. In priniple one can find for every 
SM extension designed to explain the observed neutrino data a supersymmetric version. 
The most popular among them are the various seesaw models [21, 22, 23, 24, 25, 26] leading 
to an effective dim-5 operator generating Majorana masses for the neutrinos [27, 28, 29]. 
Moreover, there is an intrinsic supersymmetric mechanism to explain neutrino data, namely 
the breaking of R-parity via Lepton number violating interactions [30, 31, 32, 33]. Another 
reason is the strong CP problem which is also still present in the MSSM. One way to solve 
this is to postulate an additional U(l) group as proposed by Peccei and Quinn [34]. In 
SUSY, this would lead to the presence of a pseudo scalar Axion, the scalar Saxion and the 
fermionic Axino [35]. Another intrinsic problem of the MSSM is the so called /U problem: 
above the breaking scale, the fj, parameter in the superpotential of the MSSM is the only 
dimensionful parameter in the model. The natural scale of that parameter would therefore 
be either of the GUT scale or exactly if it is forbidden by symmetry. However, we know 
that it must be of order of the electroweak scale to explain precision data [36] and to fullfill 
the existing bounds from collider searches [37] . One idea to solve the fi problem is to create 
an effective \i term just after SUSY breaking like in the NMSSM [38]. 

Of course, this list can be still extended, and over the years, people have come up 
with many different ways of modifying and extending the MSSM. However, the way from 
the first idea about a SUSY model to numerical results is normally long and exhaustive: 
all analytical expressions for masses, vertices and renormalization group equations have to 
derived. Code has to be generated to calculate the numerical values for the masses and, if 
the results are to be reasonable, to add loop corrections. Before the existing software tools 
to calculate widths or cross sections as well as perform Monte Carlo studied can be used, 
the new model has to be implemented. That demands not only a good knowledge about 
the different programs but is also a very time consuming task. 

In this paper, we present a framework which covers all of the aforementioned steps in 
an automatized way and which goes also further. Based on SARAH [39, 40, 41], a tool chain 
can be created for a large variety of SUSY models which covers the spectrum calculation 
with SPheno [42, 43], the calculation of cross sections with CalcHep [44, 45] and the 
determination of the dark matter relic density using micrOMEGAs [46], the analysis of 
constraints coming from the Higgs searches via HiggsBounds [47, 48] and Monte Carlo 
simulations of collider observables with WH IZARD [49, 50]. The basic idea is that the user 
can implement the model in SARAH in an intuitive and fast way. Afterwards, SARAH 
generates Fortran source code for SPheno as well as model files for CalcHep, O'Mega 
and WH IZARD. Since the implementation of the model in all of the programs is based 
on the one implementation in SARAH, the conventions are the same for all programs, 
greatly simplifying the transfer of information between the different tools. For example, 
the spectrum file of SPheno can be directly used by micrOMEGAs to calculate the relic 
density. In addition, the new SPheno modules include the possibility to write input file 
for HiggsBounds. 
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We start with an introduction to SARAH in sec. 2, demonstrating how models can be im- 
plemented and what information can be derived by SARAH. In sec. 3, we discuss the details 
of the Fortran output of SARAH for SPheno, and in sec. 4, we show how HiggsBounds 
can be applied to the model. Afterwards, the link to CalcHep and micrOMEGAs is dis- 
cussed in sec. 5, and the output for O'Mega and WHIZARD is presented in sec. 6. Finally, 
we introduce the package SSP which is designed for performing parameter scans using the 
tools supported by SARAH in sec. 7 before explaining all steps necessary to combine SARAH 
and SSP together with those programs into a closed tool chain in sec. 8. 

2. Building new models with SARAH 

2.1 Overview 

SARAH is a package for Mathematica version 5.2 or higher and has been designed to 
handle every N = 1 SUSY theory with an arbitrary direct product of SU(n) and/or U(l) 
factors as gauge group. The chiral superfields can transform under arbitrary, irreducible 
representations with regard to this gauge group, and all possible renormalizable superpo- 
tential terms are supported. There are no restrictions on either the number of gauge group 
factors, the number of chiral superfields or the number of superpotential terms. Further- 
more, any number of symmetry breakings or field rotations is allowed. A schematic picture 
of the different steps performed by SARAH is shown in Fig. 1. 

2.2 Download, installation and first evaluation 

SARAH can be downloaded from 

http: //projects .hep forge .org/ sarah/ 

The package should be extracted to the application directory of Mathematica, 
$ HOME / .Mathematica /Applications/ 
on a Unix system or 

[Mathematica-Di rectory] \AddOns\Appli cat ions \ 

in a Windows environment ($HOME and [Mathematica-Directory ] should be substi- 
tuted with the home and Mathematica installation directories respectively). 

Initially, the package itself consists of three directories: the directory Package con- 
tains all Mathematica package files, while in the directory Models the definitions of the 
different models are located. The third directory LaTeX contains DT^X packages which 
are needed for the appropriate output. During execution, a fourth directory called Output 
is generated by SARAH where the results of the different calculations as well as the model 
files for the diagram calculators are stored. 

A comprehensive manual (sarah . pdf ) is included in the package archive and can also 
be found on the web page and on the arXiv [39]. In addition, a file (models .pdf) with 
information about all models delivered with the package is part of the archive. Furthermore, 
a file with a short introduction to the main commands is included (Readme . txt) as well 
as an example illustrating their use (Example . nb). 
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Figure 1: Schematic graph of the different steps performed by SARAH. The user has access to 
the calculated information shown in boxes. The ellipses show the output which can be created by 
SARAH. 



After the installation, the package is loaded in Mathematica via 
mm: << " sarah-3 . / SARAH . m" 
Subsequently, the model is initialized by 
in [2]: Start [ "Modelname" ] ; 

where Modelname is the name of the corresponding model file. As an example, the com- 
mand would read 
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in [2]: Start ["MSSM"] ; 

for the minimal supersymmetric standard model and 
m [2]: Start [ "NMSSM" , " CKM" ] ; 

for the next-to-minimal supersymmetric model in CKM basis. 
2.3 Building new models with SARAH 

As an example of how new models can be defined in SARAH, we briefly go through the 
existing MSSM model file and then show how this model can be augmented to obtain an 
implementation of the NMSSM (see [38] and references therein). 

2.3.1 Model file for the MSSM 

In the following, we showcase the different components of the MSSM implementation in 
SARAH (see [41] for a summary of our conventions). 

1. First, give a string as internal name for the model. Make sure to use only numbers 
and letters and no spaces, because this name is later on used to identify the model 



in SPheno or CalcHep. 



ModelName = "MSSM"; 
ModelNameLaTeX ="MSSM"; 







ModelNameLaTeX provides much more freedom: any symbol and spaces can be used 
and even F/T^Xsyntax is supported. 



2. The gauge sector is £7(1) x SU (2) x SU(3) and is defined by declaring the correspond- 



ing vector superfields. 



Gauge [ 


[1] 


]={B, 


u 


1] , 


hyper charge, 


gl, 


False } ; 


Gauge [ 


[2] 


]={WB, 


su 


2] , 


left, 


g2, 


True } ; 


Gauge [ 


[3] 


]={G, 


su 


3] , 


color, 


g3, 


False } ; 



First, the name of the vector superfield is given. The second entry defines the dimen- 
sion of the group, the third is the name of the gauge group and the forth gives the 
name of the corresponding gauge coupling. If the last entry is set to True, the sum 
over the group indices is expanded and the component fields are distinguished; oth- 
erwise, the sum is left implicit. In the above example, the color charges are written 
as indices, while the sum over the isospin multiplets is expanded. 

Note that SARAH automatically adds a soft-breaking gaugino mass for every vector 
superfield. 

3. The next step is to define the matter sector. That's done by the array Fields. The 
conventions are the following. First, the root of the names for the component fields 
is given (e.g. X): the derived names of the fermionic components start with F in 
front (i.e. FX), while for scalars a S is used (i.e. SX). At second position the number 
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of generations is defined and the third entry is the name of the entire superfield. 
The remaining entries are the transformation properties with respect to the different 
gauge groups. 

Using these conventions, the doublet superfields q,l,H d ,H u are added by 



Fields [[1]] = 


{{uL, dL}, 3, q, 1/6, 


2, 


3}; 


Fields[[2]] = 


{{vL, eL}, 3, 1, -1/2, 


2, 


i}; 


Fields[[3]] = 


{{HdO, Hdm}, 1, Hd, -1/2, 


2, 


l}; 


Fields [[4]] = 


{{Hup, HuO}, 1, Hu, 1/2, 


2, 


l}; 


While for the sinj 


det superfields d c ,u c , e c 






Fields [[5]] = 


{conj[dR], 3, d, 1/3, 1, 


-3} 




Fields [[6]] = 


{conj[uR], 3, u, -2/3, 1, 


-3} 




Fields [[7]] = 


{ con j [eR] , 3, e, 1,1, 


1} 





is used. 

Note that for scalars SARAH also adds the soft masses automatically. 
4. The superpotential of the MSSM is 

W = u c Y u qH u - d c Y d qHd ~ e c Y e lH d + fiH u H d 
and represented in SARAH by 



(2.1) 



SuperPotential = { {{1, Yu } , { u, q, Hu } } , { { -1 , Yd} , { d, q, Hd} } , 

{ {-1, Ye}, {e, l,Hd} }, { { 1, \ [Mu] } , {Hu, Hd} } }; 



5. There are two different sets of eigenstates: the gauge eigenstates before EWSB and 
the mass eigenstates after EWSB. The internal names are 



NameOf States= { GaugeES , EWSB}; 



6. The gauge fixing terms for the unbroken gauge groups are 

DEFINITION [GaugeES] [GaugeFixing] = 
{ {DerfVWB], -1/(2 RXi [W] ) } , 
{Der [VG] , -1/ (2 RXi [G] ) } } ; 



(2.2) 



This corresponds to 

2t,W ^9 

The gauge fixing terms are used for the calculation of the Ghost interactions in general 
i?£ gauge. These are used to write the vertices for FeynArts depending on gauge 
fixing constants, while the CalcHep output is restricted to unitary and 't Hooft 
gauge. Internal computations like those of the onedoop self-energies are performed 
in 't Hooft gauge by SARAH. 
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7. The vector bosons and gauginos rotate after EWSB as follows 



DEFINITION [EWSB] [ GaugeSector ] = 

{ { { VB , VWB [ 3 ] } , { VP , VZ } , Z Z } , 

{ {VWB [1] , VWB [2] } , {VWm, con j [VWm] } , ZW} , 

{ { fWB [1] , fWB [2] , fWB [3] } , { fWm, fWp, fWO } , ZfW) 



This encodes the common mixing of vector bosons and gauginos after EWSB 










fw-\ 


X W,2 


= z w 


w+ 


\ V,3 J 







(2.3) 

The mixing matrices can easily be parameterized in SARAH to have the standard form 



cos @w - sin Qw | z w _ 




7S 71 ^ 




1/ 



(2.4) 



The neutral components of the scalar Higgs receive vacuum expectation values (VEVs) 
Vd/v u and split into scalar and pseudo scalar components 



H 



o 



1 



U V2 
This is encoded in SARAH by 



(vd + icrd + 4>d) , 



V2 



(v u + icr u + <j) u ) 



(2.5) 



DEFINITION [EWSB] [VEVs]= 

{ {SHdO, {vd, 1/Sqrt [2] } , { sigmad, I/Sqrt [2] } , {phid, 1/Sqrt [2] } 
{SHuO, {vu, 1/Sqrt [2] } , { sigmau, I/Sqrt [2] } , {phiu, 1/Sqrt [2] } 



9. After EWSB the particles to new mass eigenstates 



DEFINITION [EWSB] [MatterSector ] = 

{ 



{SdL, SdR 


}, {Sd, ZD}} 


{SuL, SuR 


}, {Su, ZU}} 


{SeL, SeR 


}, {Se, ZE}} 


{SvL 


}, {Sv, ZV}} 


{phid, phiu 


}, {hh, ZH}} 


{sigmad, sigmau 


}, {Ah, ZA}} 


{SHdm, conj[SHup] 


}, {Hpm, ZP}} 


{fB, fWO, FHdO, FHuO 


}, {LO, ZN}} 


{{fWm, FHdm}, {fWp, 


FHup}}, {{Lm 



{ {FeL}, 
{ { FdL } , 
{ {FuL}, 



U}, {Lp,V}}}, 
{conj [FeR] } } , { {FEL, ZEL} , {FER, ZER} } } , 
{conj [FdR] } } , { {FDL, ZDL } , {FDR, ZDR } } } , 
{conj [FuR] } } , { {FUL, ZUL} , {FUR, ZUR} } } }; 
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This defines the mixings to the mass eigenstates: first, a list with gauge eigenstates 
is given, followed by the name of the new mass eigenstates and the mixing matrix. 
Hence, the first line is interpreted as 

3 3 
d~L,ia = ^2 Z ji'*dja , d,R,ia = ^ Z-^dja (2.6) 

3=1 3=1 

while the 8th line defines the mixing in the chargino sector 



(2.7) 



10. The new gauge fixing terms after EWSB are 



_ J_ \& w - + iiw _M w G~\ 2 - ^-\d^\ 2 . (2.8) 



That reads in SARAH 


DEFINITION [EWSB] [ GaugeFixing ] = 






{ {Der [VP] , 


- 1/(2 RXi[P]) }, 


{Der[VWm]+ I Mass [VWm] RXi [W] Hpm [ { 1 } ] , 


- l/(RXi[W] 


) }, 


{Der[VZ] - Mass[VZ] RXi[Z] Ah [ { 1 } ] , 


- 1/(2 RXi [ 


Z] ) }, 


{Der [VG] , 


- 1/(2 RXi [ 


G] ) } }; 


Based on this definition, A® and Hf are interpreted in 


all calculations 


as Goldstone 


bosons. 






No particles should be integrated out or deleted 






IntegrateOut= { } ; 






DeleteParticles= { } ; 






The Dirac spinors for the mass eigenstates are 


DEFINITION [EWSB] [ DiracSpinors ] = { 






Fd - > { FDL, con j [FDR] } , 






Fe -> {FEL, con j [FER] } , 






Fu -> {FUL, con j [FUR] } , 






Fv -> { FvL, 0}, 






Chi -> {L0, con j [L0] } , 






Cha -> {Lm, conj[Lp]}, 






Glu -> {fG, conj [fG] } 

}; 
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That leads to the replacements 

'-(£)--H(£r)-'-($) 

when going from four- to two-component formalism. 
2.3.2 Creating a model file for the NMSSM 

Only a few changes are necessary to turn the above MSSM model definition into a full- 
fledged implementation of the NMSSM: 

1. Add a gauge singlet superfield 

Fields [[8]] = {sR, 1, s, 0, 1, 1 } ; 



2. Change the superpotential 

SuperPotential = { {{1, Yu } , { q, Hu, u } } , { { -1 , Yd} , { q, Hd, d} } , 

{{-l,Ye}, {l,Hd,e}}, 
{ {1, \ [Lambda] } , { Hu, Hd, s } } , 
{ {1/3, \ [Kappa] }, {s,s,s}}}; 



3. Give a VEV to the scalar component of the gauge singlet 
DEFINITION [EWSB] [VEVs]= 
{ . . ., 

{SsR, {vS, 1/Sqrt [2] } , { siqmaS, I/Sqrt [2] } , {phiS, 1/Sqrt [2] } } } ; 



4. Mix the scalar part of the gauge singlet with the Higgs and the fermionic part with 
the neutralinos 

DEFINITION [EWSB] [MatterSector ] = 
{ . . ., 

{{phid, phiu, phiS}, {hh, ZH}}, 

{{siqmad, siqmau, siqmaS } , {Ah, ZA}}, 

{{fB, fWO, FHdO , FHuO,FsR}, {L0, ZN}},... } ; 



2.4 Calculations performed by SARAH 

When a model is initialized using the Start command, it is first checked for gauge anoma- 
lies and charge conservation. If any of those checks fail, a warning is printed. Afterwards, 
the calculation of the complete Lagrangian at tree- level begins, and several tree-level results 
can be obtained after it has finished. 

Masses and tadpole equations The masses and tadpole equations are derived auto- 
matically during the evaluation of a model. The user has access to both pieces of informa- 
tion through the command MassMatrix [Particle] for the mass matrix of Particle 
and TadpoleEquation [VEV] for the tadpole equation of the corresponding field. 
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Vertices SARAH can be instructed to either calculate all vertices present in the model 
or to extract only those for specific combinations of external particles. The latter task is 
performed by 

Vertex [{Particles}, Options ] ; 

(the argument of this function being a list of external particles) , while all vertices for a set 
of eigenstates can be calculated via 

MakeVertexList [Eigenstates, Options] ; 

This searches for all possible interactions present in the Lagrangian and creates lists for the 
generic subclasses of interactions. For more details about the calculation of a supersym- 
metric Lagrangian and the conventions for extracting the vertices we refer to the appendix 
of the SARAH manual [39] . 

Renormalization group equations SARAH calculates the RGEs for the parameters of 
the superpotential, the soft-breaking terms, the gauge couplings at one- and two-loop level 
and the VEVs. This is done by using the generic formulas of [51]. In addition, to handle 
the case of several, abelian gauge groups, the rules given in [52] are implemented. The 
calculation of the RGEs can be started after the initialization of a model via 

CalcRGEs [Options] ; 

Loop Corrections SARAH calculates the analytical expressions for the one-loop correc- 
tions to the one- and two-point functions (tadpoles and self energies of all particles). These 
calculations are performed in DR-scheme using the 't Hooft gauge. This is a generalization 
of the calculations for the MSSM presented in [53]. The command to start the calculation 
is 

CalcLoopCorrections [Eigenstates] ; 
2.5 Export to external programs 

SARAH can export the information derived from the model definition in a form suit- 
able for use with a number of external programs, among them FeynArts/FormCalc, 
CalcHep/CompHep, WHIZARD/O' Mega, SPheno and also plain LT^X. In the next sec- 
tions, we proceed with a detailed discussion of this functionality and demonstrate how it 
can be used in order to facilitate a comprehensive analysis of new models. 

3. Spectrum calculation with SPheno 
3.1 Introduction to SPheno 

SPheno [42, 43] is a F95 program designed for the precise calculation of the masses of 
supersymmetric particles. For this the formulas of [53] for the 1-loop masses have been 
extended to account for the flavour structures, e.g. to calculate the 1-loop corrected 6x6 
mass matrices for squarks and charged leptons and the 1-loop corrected 3x3 mass matrix 
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for sneutrinos [54, 55]. Here the complete one-loop contributions including general flavour 
mixing and general CP-phases are taken into account. In addition, the dominant two- 
loop contributions for the Higgs boson masses are included based on [56, 57, 58, 59, 60] 
Moreover, the decay rates for 2- and 3-body decays of supersymmetric particles and Higgs 
bosons are calculated. The current version of SPheno can be downloaded from 

http://projects. hep forge. org/ spheno/ 

3.2 Combining SPheno and SARAH 

SARAH is based on Mathematica and therefore it is not sensible to do exhaustive numerical 
calculations in SARAH's native environment. As opposed to that, SPheno provides fast 
numerically routines for the evaluation of the RGEs, calculating the phase space of 2- and 
3-body decays as well as Passarino Veltman integrals and much more. Since these routines 
are model independent, they can be used for all SUSY models implemented in SARAH. 
Our approach for combining all advantages of SPheno and SARAH in order to create a very 
efficient and easy way from model building to numerical results is depicted in Fig. 2: the 
model is defined in SARAH in the usual way. SARAH calculates all analytical expressions 
needed for a complete analysis of the model. This information is exported to Fortran code 
in a way suitable for inclusion in SPheno. This generates a fully functional version of 
SPheno for the new model without any need to change the source code by hand. 

The user has control over the properties of the generated SPheno version by means 
of a special input file for SARAH. First, it is possible to define the free parameters of the 
model. Those build later on the Block MINPAR in the LesHouches input file. Second, the 
boundary conditions at the GUT-, SUSY- and electroweak scales as well as at possible 
threshold scales can be set. Third, the parameters which are to be fixed by the solutions of 
the tadpoles equations can be defined. An approximate solution to the tadpoles can also 
be given, if there isn't an analytic one. 

SARAH produces replacements for all model dependent files of SPheno. These files have 
to be copied to a new subdirectory of the SPheno directory. A Makefile for compiling 
the new model afterwards as well as a template for a LesHouches input file are written by 
SARAH. 

The command to automatically calculate all necessary information like vertices and 
RGEs and generate the source code is 

in [3]: MakeSPheno [Options ] ; 

The name of the SPheno specific input file of SARAH can be given as an option. This 
offers the possibility to easily create SPheno versions for the same model with changed 
boundary conditions or another set of free parameters. 

3.3 Features of the generated SPheno version 

Calculation of the mass spectrum The SPheno version generated by SARAH calcu- 
lates the complete mass spectrum using 2-loop RGEs and 1-loop corrections to the masses, 
including the full momentum dependence of all loop integrals. In addition, for MSSM-like 
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Higgs sectors, the known two loop corrections to the Higgs masses and tadpoles can be 
included. All calculations are performed with the most general flavor structure and allow 
for the inclusion of CP phases. 

The calculation of the mass spectrum happens in an iterative way: the gauge and 
Yukawa couplings are fitted at Mz- Afterwards, a run to the GUT scale is performed and 
the boundary conditions are enforced. The parameters are than evolved down again to the 
SUSY scale and the mass spectrum is calculated. These steps are iterated until the all 
masses have converged to a given relative precision, by default 10~ 4 . 

The routines of SARAH for calculating the gauge and Yukawa couplings closely follow 
the procedures described in [42]. The values for the Yukawa couplings giving mass to 
the SM fermions and the gauge couplings are determined at the scale Mz based on the 
measured values of the quark and lepton masses, the mass of the Z-boson Mz, the Fermi 
constant Gf and the electromagnetic coupling in the Thompson limit a em (0). The 1-loop 
corrections to the mass of W- and Z-boson as well as the SUSY contributions to muon-decay 
are taken into account into the calculation. In addition, we include the complete 1-loop 
corrections to the self-energies of SM fermions [53] and re-sum the tan/3 enhanced terms 
in the calculation of the Yukawa couplings of the 6-quark and the r-lepton as described in 
[42]. The vacuum expectation values and v u are calculated with respect to the given 
value of tan (3 at Mz- 

SUSY scale input It is also possible to define the full set of free parameters (i.e. gauge 
couplings, VEVs, superpotential and soft-breaking parameters) of the model at a specific 
scale without RGE running. These parameters are afterward used to calculate the loop 
corrected mass spectrum and the decays. To use this option with SPheno, the LesHouches 
input file 1 must contain 

Block MODSEL # 

10 # Low scale input 

12 2000. # Renormalization scale 

In this example all parameters are declared to be renormalized at 2 TeV. If the no explicit 
renormalization scale is defined by the flag 12,1 TeV is used. In addition, the numerical 
values of all parameters have to be given in the LesHouches input file, e.g. 

Block GAUGEIN # 

1 0.384499E+00 # gl 

2 0.647209E+00 # g2 

3 1.121000E+00 # g3 
Block YUIN # 

1 1 8 . 57631113E-06 # Y_u(Q)~DRbar 

2 2 3 . 63064589E-03 # Y_c (Q) "DRbar 

3 3 8 . 52980570E-01 # Y_t (Q) "DRbar 



1 We extend and partly depart from the SUSY LesHouches convention [61, 62] to obtain a greater 
flexibility for the implementation of new models. 
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Calculation of decay widths The generated version of SPheno is capable of calculating 
the widths and branching ratios of all non-SM scalars and fermions. While only two body 
decays are taken in account for decaying scalars and vector bosons, the fermion widths also 
include the three body decays to all-fermion final states. In the case of Higgs fields, the 
decays into two-photon and two-gluon final states are included at leading order [63] , as are 
the gluonic QCD corrections for the Higgs decays into quarks and squarks. [64, 65]. 

3.4 Defining the properties of the generated SPheno version 

An additional file which defines the properties of the generated SPheno function is required 
by SARAH for generating SPheno output. This file must be located in the same directory 
as the other models files for SARAH. By default, it is assumed that this file is named 
SPheno. m, but it is also possible to use other file names (see sec. 3.6). This way it is 
easily possible to generate different SPheno implementations of the same model. 
The content of the SPheno specific input file for SARAH is the following: 

1. MINPAR: A list of parameters which should be read by SPheno from the block 
MINPAR in a LesHouches file. First, the number in the block is defined, afterwards 
the variable. For example: 

MINPAR = { {l,mO} , 
{2,ml2}, 
{3, TanBeta} , 
{4, SignMu} , 
{ 5 , Azero } } ; 



Later, the values of those parameters can be communicated to SPheno by using an 



input file which contains 



Block 


MINPAR # 






1 


7 . 000000E + 01 


# 


m_0 


2 


2 . 500000E + 02 


# 


M_l/2 


3 


1 . 000000E + 01 


# 


Tan (beta) 


4 . 


1 . 000000E + 00 


# 


Sign (mu) 


5 . 


. 000000E + 00 


# 


A_0 



2. EXTPAR: It is also possible to define additional parameters for the block EXTPAR of 



the LesHouches input file by 



EXTPAR = { {Nrl, 


Varl } , 


{Nr2, 


Var2 } , 


. . . }; 





For instance, in order to give three additional VEVs as input, we can use 
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EXTPAR = { { 100, vl } , 
{101, v2}, 
{102, v3}}; 



and set the values later on in the input file by 

Block EXTPAR # 

100 1.000000E-04 # v_l 

101 1.500000E-04 # v_2 
103 2.000000E-04 # v_3 



Note that there are no hard coded entries for MINPAR or EXTPAR. This makes it 
necessary to define these blocks also for models with already existing SLHA conven- 
tions. However, this also provides more freedom in varying the model and the free 
parameters. 

3. RealParameters: By default, all parameters defined in MINPAR or EXTPAR as- 
sumed to be complex, i.e. it is possible to use also the block IMMINPAR to define 
the imaginary part. However, some Fortran functions like sin can't be used with 
complex numbers, therefore is is necessary to define parameters like tan /3 explicitly 
as real, e.g. 

RealParameters = {TanBeta}; 



4. ParametersToSolveTadpoles: For each field which can obtain a VEV, SARAH 
derives the corresponding minimum condition. These equations constrain as many 
parameters as there are VEVs in the model. ParametersToSolveTadpoles de- 
fines which parameters the tadpole equations will be solved for. 

For example, to use the standard choice in the MSSM /j,, B^, the entry reads: 
ParametersToSolveTadpoles = {\[mu], B[\mu]}; 



SARAH uses the Solve command of Mathematica to solve the tadpole equations 
for the given set of parameters. If the solution is not unique because a parameter X 
appears squared, SARAH solves the equations for the absolute squared. The phase is 
then defined by the automatically generated variable SignumX, which is expected to 
be given as input. That's for instance the case of the \i parameter in the MSSM. 

The solutions for the tadpole equations are applied by SPheno during the numerical 
analysis at the SUSY as well at the electroweak scale. For that purpose the running 
values of all parameters including the VEVs are taken as input at the considered 
scale. 

5. UseGivenTapdoleSolution: In cases, in which Mathematica won't find an 
analytical solution for the tadpole equations for the given set of parameters, this 
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variable has to be set to True and an approximated solution can be given. These 
solutions are defined by 

• SubSolutionsTadpolesTree: For the solution at tree level 

SubSolutionsTadpolesTree = {xl -> soil, x2 -> sol2,...}; 

Here, xl, x2 are the names of the parameters which are fixed by the tadpole 
equations and soli, sol2 are the approximated expressions for them. 

• SubSolutionsTadpolesLoop: The solutions of the one loop corrected tad- 
pole equations. The one loop corrections to the different VEVs have to be named 
TadlLoop [i] . 

6. RenormalizationScaleFirstGuess: For the first run of the RGEs, before any 
mass has been calculated by SPheno, the squared renormalization scale can be de- 
fined by this entry. For example, for a mSugra scenario the common choice is 

RenormalizationScaleFirstGuess = m0~2 + 4 ml2~2; 



This affects the running only if the SUSY scale is not fixed and SPA conventions are 
disabled in the LesHouches input file. 

7. RenormalizationScale: For all further runs, another renormalization scale can 
be given which is a function of the calculated masses, e.g. 

RenormalizationScale = MSu [ 1 ] *MSu [ 6 ] ; 



8. Two loop contributions to the Higgs masses: if the Higgs sector of the model is the 
same as for the MSSM, the original SPheno routines for calculating the two- loop 
tadpole equations and two-loop self energies to the the scalar and pseudo scalar Higgs 
can be activated by setting 

UseHiggs2LoopMSSM = True; 



9. Condition for the GUT scale: to set a condition for a dynamically adjusted GUT 
scale, use 

ConditionGUTscale = l.h.s == r.h.s; 



A common choice would be the unification point of g\ and </2- in that case the 
condition reads 

ConditionGUTscale = gl == g2 ; 



Note, that the value of the left-hand side must be smaller at scales below the GUT 
scale as the right hand side. 
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10. Boundary Condition: It is possible to define boundary conditions at three different 
scales: 



• Electroweak scale: BoundaryEWSBScale 

• SUSY scale: BoundarySUSYScale 

• GUT scale: BoundaryHighScale 

In addition, if thresholds are involved, boundary conditions can be set at the thresh- 
old scale, see sec. 3.7. It is also possible to use a low scale input without any 
RGE running. In that case special boundary conditions can be defined by the array 
BoundaryLowScale Input. 

All boundaries are defined by a two dimensional array. The first entry is the name 
of the parameter, the second entry is the used condition at the considered scale. The 
condition can be . . . 

• ... an input parameter from MINPAR or EXTPAR, e.g. 

{MassB, m!2}; 



• ... a block in the SLHA input file, e.g. 
{ Yv, LHInput [Yv] } ; 



• ... a function of different parameters, e.g. 
{ TYd, Azero*Yd}; 



• ... a diagonal matrix, e.g. 

{md2, DIAGONAL m(T2}; 



• . . . matrix multiplications or the inverse of a matrix, e.g. 



{X, MatMul2 [A, InverseMatrix [B] , FortranFalse ] } ; 







For the matrix multiplication MatMul2 has to be used. The third argument con- 
trols whether if only diagonal elements (FortranTrue) should be considered 
or not ( FortranFalse). 

• ... a self defined function 



{X, Func [A, B, C] } ; 







It is also possible to use some self defined function. The Fortran code of that 
function has to included in the array Self Def inedFunctions in SPheno . m. 
Later on it will be written to Model_Data . f 90. Note, that the standard func- 
tions needed for GMSB are already included [66]: 
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- fGMSB [X] : 

m = ^ (mi+«) - + I^j^j) + 

^(tp-^-^-i.j+lu.P-i.)) (3.1) 

- gGMSB [X] : 

g(x) = ln(l + x) + ln(l - x) (3.2) 



Boundary conditions can be overwritten by assigning a value to a parameter in the 
LesHouches input file. For example, the Higgs soft breaking masses at the GUT scale 
can be forced to have specific values instead of by declaring 



Block MSOFTIN 


# 




21 10000.000 


# 


mHd2 


22 20000.00 


# 


mHu2 



in the SLHA file. 



Several sets of boundary conditions In order to implement different versions of 
a single model which differ only by the used boundary conditions, BoundaryEWSBScale, 
BoundarySUSYScale, BoundaryHighScale can be also a nested list, e.g. 



BoundarySUSYScale = Table 


[ { } , { 2 } ] ; 




BoundaryGUTScale = Tablet 


{ }, {2}] ; 




BoundarySUSYScale [ [1] ] = 


{ { KappaNMSSM, 


Kappalnput } , 




{LambdaNMSSM, 


Lambdalnput } } ; 


BoundaryGUTScale [ [1] ] 


{ }; 




BoundarySUSYScale [ [2] ] = 


{ }; 




BoundaryGUTScale [ [2] ] 


{ {KappaNMSSM, 


Kappalnput } , 




{LambdaNMSSM, 


Lambdalnput } } ; 



In the first case, the input values for A and k are taken at the SUSY scale, in the 
second one at the GUT scale. To communicate to SPheno which set of boundary 
conditions should be used for a run, flag 2 in MODSEL is used: 



Block MODSEL # 
2 X # This uses the X. set of boundary conditions. 

The default value is 1. 

11. Lists for calculating decay widths: 

• ListDecayParticles: List of particles for which two-body decays will be 
calculated. This can be a list of particles using the names inside SARAH, e.g. 
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ListDecayParticles = { Sd, Su, Se, hh, Ah, Hpm, Chi } ; 

or just Automatic. If Automatic is used, the widths of all particles not 
defined as standard model particles as well as the top width are calculated by 
SPheno. 



• ListDecayParticles3B; Three body decays of fermions. This can be a list 
with the names of the particles and the corresponding files names, e.g. 



ListDecayParticles3B = 


{ {Chi, 


"Neutralino. f 90" } , 




{Cha, 


"Chargino . f 90" } , 




{Glu, 


"Gluino. f 90" } } ; 



or just Automatic. If Automatic is used, the widths of all fermions not 
defined as standard model particles are calculated. The auto generated file 
names are ParticleName . f 90. 



12. Ordering of mass eigenstates: normally, all particles of one kind are ordered in 
SPheno by their mass. However, it might be desirable to override this behavior 
and instead define another ordering scheme. For example, consider several massless 
CP odd particles at tree level exist which can be assigned to a Goldstone boson. For 
this purpose, a condition can be defined by using 

ConditionForMassOrdering = { {Particle, Condition}, ... }; 



The condition has to be Fortran source code and is added to the corresponding 
routine. For instance, a condition for the NMSSM would read 



ConditionForMassOrdering= 


= { 


{Ah, 




"If (Abs (ZA (1, 3) ) . gt . Abs (ZA (1, 2) ) ) Then \n 


MAh2temp = MAh2 \n 




ZAtemp = ZA \n 




ZA (1, : ) = ZAtemp (2, : ) 


\n 


ZA (2, : ) = ZAtemp (1, : ) 


\n 


MAh2(l) = MAh2temp(2) 


\n 


MAh2(2) = MAh2temp(l) 


\n 


End If \n \n" } 

}; 





This example checks whether two massless pseudo scalars are present in the spectrum 
and, if this is the case, uses the not singlet-like particle as Goldstone boson. 

3.5 Example: The input file for the MSSM 

To generate the SPheno output for the MSSM, SPheno. m should include the following 
information: 
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1. We want to have mSugra like boundary conditions. Therefore, we chose the minimal 
set of parameters defining the model as mo, M 1 / 2 , A)> sign^i and tan/?. These will 
later be read from the MINPAR block of a LesHouches input file. 

MINPAR={ {l,mO}, 
{2,ml2}, 
{3, TanBeta} , 
{ 4 , SignumMu } , 
{ 5 , Azero } } ; 



2. TanBeta has to be declared as a real parameter 
RealParameters = {TanBeta}; 



3. As usual in the MSSM, the tadpole equations should be solved with respect to ji and 



ParametersToSolveTadpoles = { \ [Mu] , B [ \ [Mu] ] } ; 



4. To study models with a dynamically adjusted SUSY scale, the expressions for the 
definition of the SUSY scale can be given. The first expression is used only before 
the mass spectrum has be calculated the first time. Note, that these definitions can 
easily disabled in the LesHouches input file by flag MODSEL 12 and a fixed scale can 
be used. Also, when SPS conventions are switched on in the LesHouches input file 
by SPheno Input 2, a fixed scale of 1 TeV is used. 

RenormalizationScaleFirstGuess = m0~2 + 4 ml2~2; 
RenormalizationScale = MSu [ 1 ] *MSu [ 6 ] ; 



5. The GUT scale is the unification point of g\ and g 2 

ConditionGUTscale = gl == g2 ; 



6. As said, we want to use mSugra like boundary conditions. These are straightforward 
define by 



BoundaryHighScale= { 



{ T [Ye ] 


, Azero*Ye } , 




{ T [Yd] 


, Azero*Yd} , 




{T[Yu] 


, Azero*Yu } , 




{mq2, 


DIAGONAL mO 


~2}, 


{ml2, 


DIAGONAL mO 


"2}, 


{md2, 


DIAGONAL mO 


"2}, 


{mu2 , 


DIAGONAL mO 


"2}, 


{me 2 , 


DIAGONAL mO 


"2}, 
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{mHd2, m(T2}, 

{mHu2 , mO " 2 } , 

{MassB, ml2}, 

{MassWB, ml2}, 

{MassG, ml2} 

}; 



7. It is also possible to use the generated SPheno version with a low scale input. This is 
enabled by setting MODSEL 1 to 0. In that case, input values for all free parameters 
of the model are expected. However, also in this case, we can define a set of boundary 
conditions, e.g. for dynamically calculating the SUSY VEVs 



BoundaryLowScaleInput= { 

{vd, Sqrt [2 mz2/ (gl~2 + g2~2) ] *Sin [ ArcTan [ TanBeta ; 

(vu, Sqrt [2 mz2/ (gl~2+g2~2) ] *Cos [ArcTan [ TanBeta ; 
}; 



8. Finally, we define that the code for the calculation of the two and three body decays 
is generated for all SUSY particles. That's done by using the flag Automatic. 

ListDecayParticles = Automatic; 
ListDecayParticles3B = Automatic; 



3.6 Generating the output 

After an input file with all necessary information has been created, the generation of the 
source code for SPheno is triggered by 

in [5]: MakeSPheno [Options ] 

The different options are: 

• Eigenstates->Name of Eigenstates . If not specified, the last set of eigen- 
states is used which corresponds to the last entry of NameOf States (see sec. 2.3.1). 

• ReadLists->True can be used if all vertices and RGEs have already been calcu- 
lated for the model and the former results should be used to save time. 

• InputFile. The name of the SPheno input file. If not defined, SPheno . m is used. 
The generated source code is located in 

[SARAH Directory] /Output/ [Model] / [Eigenstates] /SPheno/ 

It is sufficient to copy all files from this directory to a sub-directory NAME of SPheno 
version 3.1.4 or later and compile it using in the SPheno root directory 



> make Model=NAME 



NAME has to be the same as the name of the model defined in SARAH (see sec. 2.3.1). 
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3.7 Including Thresholds 

Using SARAH it is possible to include thresholds in the RGE running performed by SPheno. 
3.7.1 Thresholds without gauge symmetry breaking 

If all scales have the same gauge structure, it is possible for SARAH to derive the RGEs for 
all scales from the RGEs for the highest scale by performing the following steps: 

• For those fields which should be integrated out during the run, variables n gen {^i) are 
introduced internally which define the number of generation of the heavy field <I> j . All 
gauge group constants like the Dynkin index summed over chiral superfields, S(R), 
are expressed as function of n gen (&i). These n gen {^i) are dynamically adjusted, when 
the energy scale crosses a threshold. 

• When crossing a threshold, the couplings involving heavy fields are set to zero. For 
example, the Yukawa type coupling of the form Y %3 <&i(j>jH involves three generations 
of the heavy field <F At the threshold of <E>fc, the fc-th row of Y is set to zero. That 
happens similarly for all other superpotential and soft-breaking parameters. 

• The masses of scalar and fermionic components of the heavy superfields are assumed 
to be identical, i.e. the soft SUSY breaking terms are assumed to be negligible. These 
masses are given by a bilinear superpotential term. 

In order to include thresholds without gauge symmetry breaking, the following steps 
have to be performed: 

1. The heavy fields must be marked for deletion in the SARAH model definition: 
DeleteFields = {...}; 



This ensures, that the decays, loop corrections, etc. at the SUSY scale calculated by 
SPheno receive no contributions from the heavy fields 

2. The thresholds have to be defined in SPheno. m : 

Thresholds = {{Scalel, {HeavyFieldsl } } , 
{Scale2, ... } } ; 



For all scales an entry in the array Thresholds has to be added. Each entry defines 
the threshold scale and at second position a list of the heavy superfields which can 
be restricted to specific generations. 

It is possible to define boundary conditions at each threshold scale for running up and 
down separately: 

BoundaryConditionsUp [ [x] ] = { . ..}; 
BoundaryConditionsDown [ [x] ] = { . . . } ; 
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Threshold corrections Using 2-loop RGEs requires 1-loop boundary condition. There- 
fore, at each threshold scale the one loop threshold corrections to gauge couplings and 
gaugino masses are calculated. The general expressions are [14] 

^^leb^'KCO ' (3 ' 4) 

T^if) is the Dynkin index of a field transforming as representation r with respect to the 
gauge group belonging to the gauge coupling M is the mass of this particle and Mt is 
the threshold scale. 

Example As an example, a version of SPheno implementing the seesaw type II and 
type III models can be generated by adding the following entries to Spheno .m 

1. Seesaw II: 

Thresholds= { 

{Abs [MTMIN] , {s, sb, t, tb, z, zb} } 

}; 



2. Seesaw III: 



Thresholds= { 


























{Abs [MWM3IN 


1,1] 


] , {Hx3 


:r 


, Hxb3 


:n 


,Hg3 


[1] 


, Hb3 


:i] 


, Hw3 


:i] ] 


}, 


{Abs [MWM3IN 


2,2] 


] , {Hx3 




, Hxb3 




,Hg3 


:2] 


, Hb3 


:2] 


, Hw3 


12] ] 


}, 


{Abs [MWM3IN 

}; 


3,3] 


] , {Hx3 


:3: 


, Hxb3 


:3] 


,Hg3 


[3] 


, Hb3 


:3] 


, Hw3 


13] ] 


} 



For a more comprehensive discussion of the model files for the Seesaw I, see appendix A.l. 
3.7.2 With gauge symmetry breaking 

If the gauge structure at the different scales are different, each set of RGEs is calculated 
separately and this information is then combined into one consistent version of SPheno 
which includes routines for calculating finite shifts in the gauge couplings, gaugino and 
scalar mass parameters. As an example, the implementation of a left-right supersymmetric 
model with one symmetry breaking scales is shown in sec. A. 2. 

In order to implement such a model, the following steps are necessary: 

1. For each regime, a separate model file for SARAH has to be created. These model file 
have to be saved in the subdirectories Regime- 1, Regime-2, ... of 

[SARAH Directory] /Models/ [Model] / 



beginning with the highest scale. 
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2. The SPheno input file for the higher scales must provide the following information: 

• IntermediateScale = True 

• RegimeNr = X 

• A list of the heavy fields, which should be integrated out, the gauge sector below 
the threshold as well as the corresponding quantum numbers of the fields which 
are to be integrated out. This is needed to calculate the finite shifts at the 
threshold scale, for instance 



HeavyFields = 


{Field_ 


1, Field_ 


2, ..}; 




NextGauge = 


{U[l] , 


SU[2] , 


SU[2] , 


SU[3] }; 


NextQN = { 










{Field, 


.1, o, 


2, 


1, 


1}, 


{Field. 


2, 1/3, 


1, 


2, 


4}, 


}; 











Care has to be taken to use the same ordering of the SU(2) gauge groups as in 
the orignal definition. 



3. All necessary information for combining the regimes to one SPheno is given in 
SPheno . m of the lowest scale. 

• IntermediateScale = False 

• RegimeNr = X 

• The threshold scales: ThresholdScales = ... 

• The boundary conditions for running up and down at each threshold scale: 

BoundaryConditionsUp [ [x] ] = {...}; 
BoundaryConditionsDown [ [x] ] = {...}; 

In the boundary conditions indexl, index2, . . . can be used for defining sums 
over indices. 

• The usual information for SPheno, defined in the sec. 3.4. 

When starting the SPheno output of the lowest scale, all other scaler are evolved 
automatically. Note that, in order to calculate the RGEs of the different regimes, SARAH 
starts one additional Mathematica kernel. For passing the information between the 
different Mathematica kernels a directory Meta in the model directory is created by 
SARAH. The screen output of Mathematica during the evaluation of the higher regimes 
is written to that directory (Output-Regime-X . m), allowing the user to supervise the 
progress and see potential error messages. The necessary information of each regime for 
writing the combined source code for SPheno at the end is saved by SARAH in the files 
Regime-X . m. 
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3.8 Low energy SPheno version 

It is also possible to create a SPheno version with much less features which only accepts 
low energy input. That means, the RGEs are not written out and also the fit to the 
electroweak data is not performed in the numerical evaluation of one point. It just solves 
the tadpole equations, calculates the tree- and one-loop masses as well as the decay widths 
and branching ratios. The advantage of such a SPheno version is that it works with a 
larger set of models, e.g. also non-SUSY models or other models not supported by a full 
evaluation as explained in sec. 3.10. To get a SPheno version without RGE evolution, 
insert 

OnlyLowEnergySPheno = True; 

in SPheno . m. The remaining information needed by SARAH is only a small subset of the 
settings discussed above and consists of 

• MINPAR 

• ParametersToSolveTadpoles 

• BoundaryLowScalelnput 

• Li stDecayP articles and ListDecayParticles3B. Note that the Automatic 
statement for automatically deriving the decays of all non-SM particles does not work 
in this case as SARAH doesn't differ between SUSY or Non-SUSY particle in order to 
make the output as generic as possible. Therefore, the lists of the decaying particles 
have to be supplied manually. 

3.9 Differences to SPheno3 .1.4 

A few things are handled in a slightly different way in a SPheno MSSM module created 
by SARAH in comparison to SPheno 3.1.4. Also, some things are not yet implemented. 
We give an overview in Tab. 1. 

3.10 Supported models and known issues 

While SARAH can create valid SPheno code for many different models, there are some 
requirements on the model and some minor restrictions on the functionality of the resulting 
SPheno module. At the moment, those are 

• Fit to low energy data: in order to perform a fit to low energy data (e.g. for 
fermion masses, rnz, Gf and a em ) as starting point of the RGE evaluation, the 
following parameters must be present in the model: Yukawa couplings for lepton 
and quarks, two Higgs VEVs and, of course, the three SM gauge couplings and the 
SM particle content. However, it is still possible to use at least some features of 
the SPheno output of SARAH by manually supplying model parameters for SPheno. 
In that way, the RGE evaluation and the fit the electroweak data is skipped, but 
the one-loop corrected masses as well as the decay widths and branching ratios are 
calculated. 
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SPheno3 .1.4 


SPheno by SARAH 


Loop corrections in 
Higgs sector 


Complete one-loop correc- 
tions and dominant two-loop 
corrections. 


One-loop corrections are cal- 
culated by SARAH. Dominant 
two-loop corrections imple- 
mented in SPheno can be 
linked. 


Tadpole equations 


Solved at SUSY scale. Solu- 
tions are evaluated at the EW 
due to RGE running 


Solved at SUSY and EW 
scale. 


Three body decays 


Three body decays of fermions 
and stop 


Three body decays of fermions 


Loop induced 
neutralino decays 


One-loop decay into photon 
and neutralino included 


No loop induced decays in- 
cluded 


Loop induced 
gluino decays 


One-loop decay into gluon 
and neutralino included 


No loop induced decays in- 
cluded 


e + e~ collisions 


Calculates cross sections 


Cross section calculation not 
included 



Table 1: Different handling and implementation status in SPheno 3.1.4 vs. the SPheno modules 
created by SARAH. 



• Flavor decomposition: with SARAH it is possible to assign a unique name to 
each generation of a particular field and this way treat the individual generations as 
independent fields. That is not yet supported in the SPheno output. Furthermore, 
mixing matrices generated with the option NoFlavorMixing can not yet be handled 
by the numerical code. 

4. Checking Higgs constraints with HiggsBounds 

HiggsBounds [47, 48] is a tool to test the neutral and charged Higgs sectors against 
the current exclusion bounds from the Higgs searches at the LEP, Tevatron and LHC ex- 
periments. The required input consists of the masses, width and branching ratios of the 
Higgs fields. In addition, it is either possible to provide full information about produc- 
tion cross sections in e + e~ and pp collisions, or to work with a set of effective couplings. 
HiggsBounds can be downloaded from 

http://projects. hep forge. org/higgsbounds 

Although HiggsBounds supports the LesHouches interface, this functionality is re- 
stricted so far to at most 5 neutral Higgs fields, and therefore, we don't use it. Instead, 
SPheno modules generated by SARAH can create all necessary input files needed for a run 
of HiggsBounds with effective couplings (option whichinput=ef f C). To write these 
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files, the flag 7 5 in the block SPhenoInput in the LesHouches input file has to be set to 
1. 

Block SPhenoInput # 

75 1 # Write files for HiggsBounds 

Unfortunately, we can not provide all information which can be used by HiggsBounds to 
check the constraints. In particular, the effective couplings H — > and H —> ggZ are not 
yet calculated by SPheno and therefore set to zero in the output. In addition, as already 
mentioned, the SPheno version created by SARAH does not support calculating the e + e~ 
cross sections. For this reason, also the LEP production cross section of charged Higgs 
fields is not available for SPheno and it sets this value also to 0. However, it is of course 
possible to calculate this cross section as well as all other cross sections needed for the 
options (whichinput=hadr or whichinput=part) of HiggsBounds using CalcHep 
or WHIZARD with the corresponding model files created by SARAH (cf. sec. 5 and 6). 
The following files are written by SPheno 

• MH_GammaTot.dat: 

Masses and widths of all neutral Higgs fields 

• MHplus_GammaTot.dat: 

Masses and widths of all charged Higgs fields 

• BR_H_NP . dat: 

Branching ratios of neutral Higgs fields into invisible and other neutral Higgs fields. 

• BR_Hplus . dat: 

Branching ratios of charged Higgs fields into cs, cb and tv final states 

• BR_t . dat: 

Top quark branching rations into a bottom quark plus either a W boson or a charged 
Higgs 

• ef f C . dat: Effective couplings among neutral Higgs fields and of neutral Higgses to 
ss, cc, bb, ti, up,, tt , 77, gg, jZ, ggZ 

• LEP_HpHm_CS_ratios.dat: 

LEP cross sections for the production of charged Higgs bosons (set to zero, see above) 

HiggsBounds can be run on the generated files by invoking it as 

> ./HiggsBounds LandH effC [NN] [NC] [SPheno Directory] 

where [NN] has to be replaced by the number of neutral Higgses and [NC] by the number 
of the charged ones. Please consult the HiggsBounds manual for further information. The 
results of the check are written to the file HiggsBounds_results . dat which is located 
in the same directory as the input file, i.e. in our case in the SPheno root directory. 
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5. Calculation of cross sections, widths and relic density using CalcHep 
and micrOMEGAs 

CalcHep [45, 44] is a package for the calculation of Feynman amplitudes and their inte- 
gration over multi-particle phase space. It is designed to provide a direct transition from 
the Lagrangian to the cross sections and distributions. The CalcHep homepage is located 
at 

http :/ /theory . sinp . msu . ru/ ~pukhov/calchep . html 
5.1 Using CalcHep with SARAH and SPheno 

The CalcHep model files produced by SARAH support both the Feynman and unitarity 
gauges. Furthermore, SARAH can split interactions between four colored particles as re- 
quired by CalcHep/CompHep. Models with CP violation are also supported. The model 
files for CalcHep/CompHep are created by 

in [5]: MakeCHep [Options ] ; 

Options exist for specifying the gauge (FeynmanGauge — > True /False) and for acti- 
vating the support for CP violation (CPViolation — > True). In addition, the split- 
ting of specific four-scalar interactions can be suppressed as long as they are not colored 
(NoSplitting — > list of fields) and the running of the strong coupling constant can be 
included as it is usually done in the standard CalcHep files (UseRunningCoupling — > 
True). 

Recently, it became also possible to use SLHA files with CalcHep to provide the 
numerical values of the parameters [67]. SARAH uses this option if the flag SLHAinput 
— )• True is used. Thus, the parameters calculated by SPheno can be directly passed to 
CalcHep. 

In order to calculate the same process for different points in parameter space, the so 
called blind mode of CalcHep can be used for performing calculations without the need 
to start the graphical interface. The blind mode is most conveniently put to work using 
the following steps: 

1. Start CalcHep, insert a process and choose make n_calchep. This creates an 
executable file called n_calchep in the results subdirectory of your CalcHep 
model. 

2. In order to keep a copy of the files, copy the content of results to another directory 
(make sure that the subdirectory aux is included in the copy). The spectrum file 
generated by SPheno goes into the same directory. 

3. Change to the directory containing n_calchep and start it via 
> ./n_calchep +blind 

4. Make the desired adjustments and start the numerical integration using Vegas. 
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5. When n_calchep exits, a line of the form [{]][{{[]] is returned which mimics 
the keyboard commands you did 

6. Now, the same calculation can be repeated without the CalcHep frontend using the 
command 

> ./n_calchep -blind "[{]][{{[]]" 

This way, a number n.calchep instances for several processes (e.g. sparticle or Higgs 
production at LEP, Tevatron or LHC) can be created and called non- interactively from a 
script, thus facilitating automatized scans over parameter space. 

5.2 Relic density calculations with micrOMEGAs 

micrOMEGAs [46] is a well known tool for the calculation of the relic density of a dark 
matter candidate. The download is located at 

http : / /lapth . in2p3 . f r /mi cr omegas/ 

As micrOMEGAs uses CalcHep for the calculation of (co-) annihilation cross sections, a 
model file for CalcHep must be generated first before micrOMEGAs can be used. 

SARAH writes two files for micrOMEGAs which can serve as so-called main files, i.e. 
they can be compiled with micrOMEGAs and executed to perform calculations. While 
CalcOmega . cpp calculates only Q,h 2 and writes the result to the file omg.out. 
CalcOmega_with_DDetection . cpp computes also direct detection signals. As the SLHA+ 
import functionality of CalcHep can also be used with micrOMEGAs, it is sufficient to 
simply copy the spectrum file written by SPheno to the directory of micrOMEGAs and 
start the calculation. 

6. Monte Carlo studies with WH IZARD 
6.1 Introduction 

WHI ZARD [50] is a fast tree-level Monte Carlo generator for parton level events. A particular 
strength of the code is the efficient generation of unweighted events for high multiplicity 
final states (simulations with 8 final state particles have been performed successfully) using 
exact matrix elements. This makes it particularly useful for the study of supersymmetric 
models which generically feature complicated multiparticle final states arising from long 
decay chains. 

Behind the scenes, WHI ZARD builds on the optimizing matrix element generator ' Mega 
[49]. In order to deliver the fast tree level matrix elements required by WHI ZARD, 0' Mega 
builds up matrix elements as directed acyclical graphs of one-particle off shell wave func- 
tions (1POW, Green's functions with all but one legs amputated). These graphs are then 
transformed into highly optimized FORTRAN 90 code which is called by WHIZARD to 
calculate helicity matrix elements by recursively fusing lPOWs. This algorithm is guaran- 
teed to avoid any redundancies arising from the repeated evaluation of subdiagrams from 
the start and can be shown to grow only exponentially in complexity with the number of 
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external legs (as opposed of the factorial growth of any traditional Feynman diagram based 
approach) . 

For the treatment of color, WHIZARD and O'Mega leverage the color flow decomposi- 
tion [68]. The amplitudes are decomposed into all possible different color flows, for which 
the amplitudes are then calculated using 0' Megas recursive algorithm and finally combined 
into the squared and color summed matrix elements. In addition to the fast calculation of 
the color trace, this algorithm also provides the color connection information later required 
by fragmentation and hadronization algorithms. 

The integration and event generation in WHIZARD is performed using the adaptive 
multichannel Monte Carlo code VAMP [69]. For each process, WHIZARD dynamically deter- 
mines a set of suitable phase space maps. VAMP assigns a grid to each map, all which are 
then linearly combined to form the phase space parameterization. During the integration 
process, VAMP adapts both the grids and the weights of the different channels. After in- 
tegration, the optimized grids are used to facilitate the efficient generation of unweighted 
events. 

Recently, WHIZARD has received a major upgrade, going from the (now legacy) 1.x 
version branch (currently 1.97) to the modern 2.x branch (currently 2.0.5). While the 
interface between SARAH and WHIZARD supports both branches, the new version features 
both many new physics features (including factorized matrix elements incorporating full 
spin and color correlations) and technical improvements, so using the newer version is 
highly advised. The WHIZARD package can be downloaded from 

http : / / www . hepf orge .org/ archive/whizard/ 
6.2 The interface between SARAH and WHIZARD 

The interface between SARAH and WHIZARD shares significant parts of its code with the 
interface between FeynRules [70], with a thin layer on top to interface with SARAH. As 
such, the feature set is identical with that of the FeynRules interface and, although we 
will now give a short overview, we point the reader to above reference for more information. 

At the moment, the interface can handle spin 0, g an d 1 fields, with support for spin 
| being being planned for a future revision. The supported interactions are a subset of the 
operators currently available in 0' Mega. Specifically, up to very few exceptions, nearly all 
dimension 3 and 4 interactions are available, together with a couple of higher dimension 
operators. In particular, this list is sufficient to handle all interactions which are generated 
in a typical application of SARAH to a supersymmetric model without higher-dimensional 
terms in the Lagrangian. 

The list of supported color structures is more complicated, but in practice all color 
tensors arising from QCD gauge interactions which lead to supported Lorentz structures 
are supported. Again, this covers all interactions which can be generated by SARAH. A full 
list of the supported color structures can be found in [71]. 

The 2.x branch of WHIZARD supports running the strong coupling, and the interface 
supports this functionality. In order to use it, the model must define a parameter called 
aS which is treated by the interface as as at the Z pole. In addition, the QCD coupling 
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WOWhizardVersion 


Whizard versions supported 


"1 . 92" 


1.92 


"1 . 93" 


1.93 - 1.95 


"1 . 96" 


1.96+ 


"2 .0" 


2.0 - 2.0.2 


"2.0.3" (default) 


2.0.3+ 



Table 2: Currently available version choices when generating WHIZARD model files, together with 
the respective WHIZARD versions supported by them. 

constant gs must be derived from as and have the description "Strong-Coupling". 
Once these criteria are met, the generated model will automatically evolve all vertex factors 
which depend on either gs or as if the running coupling is activated in WHIZARD. For 
example, all model files included in SARAH support the running by deriving all vertices as 
function of gs and including the following definitions into the parameter list: 



{g3, { 


• • • r 






Description -> "Strong-Coupling", 






DependenceNum -> Sgrt [Alphas 4 Pi] 


} }, 


{Alphas, { 


. . . , 






OutputName-> aS } } , 





6.3 Generating model files and using them with WHIZARD 
Calling the interface 

In order to generate model files for WHIZARD, the Model Output command is used 
in [3j: ModelOutput [Eigenstates, WriteWHIZARD->True] ; 

This writes the model file using the default options. An equivalent syntax which allows for 
passing options to the interface is 

in [3]: ModelOutput [Eigenstates ] ; 
in [4]: MakeWHIZARD [Options] ; 

The different options accepted by MakeWHIZARD are: 

1. Exclude, Values: list of generic vertex types, Default: {SSSS} 

Prevents vertices matching the generic types from being generated in order to speed 
up the program and reduce the complexity of the generated model. 

2. WOModelName, Values: string, Default: predefined model name 

Gives the possibility to change the model name. If output for WHIZARD 1.x is gen- 
erated, the name should start with fr_ in order for the model to be picked up 
automatically by the WHIZARD build system. 
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3. MaximalCouplingsPerFile, Values: Number, Default: 500 

Defines the maximal number of couplings written to one file. Adjusting this eases 
the workload on the compiler when compiling the model. 

4. Version, Values: String, Default: latest version 

Defines the version of WHIZARD for which the model file is generated. A list of all 
values for this setting currently valid can be found in Tab. 2. In addition, you can get 
a list via ?WO "WhizardVersion in Mathematica after SARAH has been loaded. 

5. ReadLists, Values: True or False, Default: False 

This setting controls whether the cached results from a previous calculation should 
be used. 

While generating the model files, the interface will print status messages and infor- 
mation about potential incompatibilities of the model with WHIZARD to the screen. It is 
highly advised to read through this information carefully. 

Using the generated model files with WHIZARD 

After the interface has completed, the generated files can be found in the WHIZARD_Omega 
subdirectory of SARAH's output directory. In order to use the model with WHIZARD 2.x, 
the generated code must be compiled and installed. For most applications, this is done by 
simply issuing (inside the output directory) 

> . / configure 

> make 

> make install 

By default, the third command installs the compiled model into . whi zard in current user's 
home directory where it is automatically picked up by WHIZARD. Alternative installation 
paths can be specified using the --prefix option to WHIZARD. 

> ./configure --pref ix=/path/to/installation/pref ix 

If the files are installed into the WHIZARD installation prefix, the program will also pick 
them up automatically, while WHIZARD's --localpref ix option must be used to com- 
municate any other choice to WHIZARD. In case WHIZARD is not available in the binary 
search path, the WO_CONFIG environment variable can be used to point configure to 
the binaries 

> ./configure WO_CONFIG=/path/to/whizard/binaries 

More information on the available options and their syntax can be obtained with the 
--help option. 

In the case of WHIZARD 1.x output, the generated files must be patched into the 
WHIZARD source tree. To this end, the interface creates a script called inject. In most 
cases, it is sufficient to simply call the script as 

> ./inject /path/to/whizard 

(from within the output directory). Issuing ./inject --help will display a list of 
options which can be used to adapt the script to more complicated usage scenarios. 
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Spectra and parameters 

In order to communicate the numerical values of the parameters calculated by SPheno to 
WHIZARD, each SPheno version generated by SARAH is capable of writing out a separate 
file providing this information as SINDARIN code which can be directly included into the 
WHIZARD input script. The user just has to make sure to use 

Block SPhenoInput # 
76 1 # WHIZARD file 

in the LesHouches input file of SPheno. The created file WHIZARD. par. [Model] can 
be included in the input script for WHIZARD by means of the include statement. For 
example, for the MSSM, the corresponding line would read 

include ( "WHIZARD . par . MSSM" ) 



7. Parameter scans with SSP 

In the previous sections, we have discussed the implementation of SUSY models into dif- 
ferent tools which can be combined to cover much of the analysis of new models. The 
corresponding work flow is depicted in Fig. 2. 

Define Model in SARAH 

I 

Derive necessary information: couplings, masses, loop 

contributions, RGEs, decay channels, . . . 

; 

Generate model dependent source code for SPheno 

I 

Compile and run SPheno 

I 

Give masses/parameters to WHIZARD, CalcHep, 
micrOMEGAs and HiggsBounds 



Figure 2: The model is defined in SARAH. Afterwards, SARAH generates all necessary files to 
implement this model in SPheno, micrOMEGAs and WHIZARD. In addition, the SPheno version 
writes a file which can be used as input for HiggsBounds. This provides an completely automatized 
way from model building to phenomenology. 

However, in order to use this tool chain it is still necessary to run all programs and 
pass the information between them. To ease this task, we have created the SSP (SARAH 
Scan and Plot) Mathematica package. This tool facilitates a quick overview of a 
new model by performing the different analysis steps automatically for many different 
points in parameter space. Furthermore, it is possible to define additional parameter space 
constraints and to use Mathematica's intrinsic functions for creating 2D parameter space 
plots. 
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7.1 Installation and running scans 

The package can be downloaded from 

http : / /projects . hepf orge .org/ sarah/ SSP . html 

After extracting the package archive to the Mathematica application directory, SSP is 
loaded via 

in[ij: <<"SSP/SSP .m" 

A scan is started by 

m [2]: Start [ " Inputf ile" ] ; 

where Inputf ile is a file containing all necessary information for a scan. 

7.2 The input files 

Each scan performed by SSP is based on the information given in two input files. In a 
first file, which has to be located in the root directory of SSP, the location of the different 
programs is given and some basic information about the different tools like the name of 
the in- and output file has to be defined. For instance, the necessary entries for SPheno 
read 

DEFAULT [SPheno] = "[SPheno Directory ] /bin/SPheno [Model ]" ; 
DEFAULT [SPhenoInputFile] = "LesHouches . in . [Model] " ; 
DEFAULT [SPhenoSpectrumFile] = " SPheno . spc . [Model ]" ; 

The second file contains all information defining a specific scan: the name of the re- 
quested settings file, the programs which should be included into the scan, input parameters 
and ranges, any constraints which should be applied and the definitions of the actual plots. 
As several different scans can be defined in one input file, each scan has first to be assigned 
an unique identifier, e.g. 

RunScans = { ScanMO , ScanM12 } ; 

In order to toggle the inclusion of the different tools, the flags IncludeWHIZARD, 
IncludeHiggsBounds, IncludeCalcHep or IncludeMicrOmegas are used 

DEFINITION [ScanMO] [IncludeHiggsBounds] = True; 
DEFINITI0N[ScanM12] [IncludeHiggsBounds] = False; 



The information on the parameters is defined by different lists for each block of the 
LesHouches input file which might look like 



DEFINITION[ScanM0M12] [MINPAR]={ 






{{1}, {Min->0,Max->1000, Steps- 


->50, 


Distribution->LINEAR} } , 


{{2}, {Min->0,Max->1000, Steps- 


->50, 


Distribution->LINEAR} } , 


{{3}, {Value->10} }, 






{{4}, {Value->1}}, 






{{5}, {Value->0}} }; 
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Other possible distributions are LOG and RANDOM. 

7.3 Features 

Fit to constraints 

SSP can be instructed to adjust one or several model parameters such that a list of con- 
straints is fulfilled. To this end, the Mathematica function NMinimize is utilized to find 
a combination of free parameters which leads to the smallest x 2 for the given constraints. 
For instance, to fix the Higgs mass between 114.9 and 115.1 GeV by a variation of Aq and 
tan/3, the necessary input is 

DEFINITION [FITHIGGS] [FitValues ] = { 
{MASS [25] , 115, . 1} }; 

DEFINITION [FITHIGGS] [FreeParameters ] = { 
{ TANBFIT, {5, 15} } , 
{AOFIT, {0,100}} }; 

DEFINITION [FITHIGGS] [MINPAR]={ 

{{1}, {Min->0 . ,Max->1000, Steps->10, Distribution->LINEAR} } , 

{{2}, {Value->500 . } } , 

{{3}, { Value -> TANBFIT } } , 

{{4}, {Value->1. } }, 

{{5}, {Value->A0FIT} } }; 

In order to tune the behavior of NMinimize, e.g. manually choosing a fit algorithm, 
options can be passed to the function 

DEFINITION [FITHIGGS] [FitOptions ] = {Method->"NelderMead" } ; 

As most of the fit algorithms cannot deal with a finite parameter range [a, b], SSP trans- 
forms to infinite boundaries via stretching 

b - a &t&n(P int ) 
P ex t = a^ 1 (b-a) (7.1) 

2 7T 

with Pi n t G [— oo,+cx)] and P ex t £ [a,b]. For more details on NMinimize and its various 
options, we refer to the Mathematica manual. 

2D parameter sampling 

One common problem is the sampling of a two dimensional parameter space, e.g. when 
checking the dark matter relic density in the (mo,M 1 / 2 ) plane. However, in many cases, 
a fixed grid or random scan might not be the best choice because certain areas should be 
sampled more precisely than others. The ContourPlot function of Mathematica was 
developed exactly for such purposes and therefore, SSP can be instructed to use it. 
The input for using ContourPlot with SSP to do a (mo,M 1 / 2 ) scan is 
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DEFINITION [M0M12] [ Count our Scan ] = 






{ DARKMATTER [ 1 ] , 






{ CONTOURSCANPARAMTER [1], 0,1500}, 






{ CONTOURSCANPARAMTER [2], 0,1500}, 






ContourPlotOptions, "DM_A0 .eps" } ; 






DEFINITION[M0M12] [MINPAR]={ 






{ {1}, 


{ Value ->CONTOURSCANPARAMTER [ 1 ] 


} }, 




{ {2}, 


{Value->CONTOURSCANPARAMTER[2] 


} }, 




{ {3}, 


{Value->10 } } , 






{ {4}, 


{Value->1 } } , 






{ {5}, 


{Min->l,Max->10 00, Steps->3 , Distribution->LOG } } 


}; 



With this setup, contour plots for the relic density in the (mo, M1/2) plane are created for 
Aq = 1, 10, 1000 GeV and tan/3 = 10, sign/x > 0. As options for the scan, the usual options 
of ContourPlot can be used. For instance, in Fig. 3, we used PrecisionGoal->"Quality " 
and varied the Plot Points between 5 and 50. The number of valid parameter points 
which have been evaluated can be found in Tab. 3. 




Figure 3: {m 0l Mi/ 2 ) plane for different values of PlotPoints. First row (from left to right): 5, 
10, 20. Second row (from left to right): 30, 40, 50. The mesh gives an insight, where Mathematica 
has decreased the distance between the interpolation points. 
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PlotPoints 


5 


10 


20 


30 


40 


50 


Calculated points 


990 


2077 


4589 


9242 


13324 


17416 



Table 3: Number of evaluated parameter space points in the (mo, Mi/ 2 )-plane for a given value of 
PlotPoints. 

8. Putting the programs to work 

8.1 Automated setup and model implementation 

In order to simplify the setup of the different programs and the generated model files, we 
provide a script to setup and install the complete environment and another script for the 
automated implementation of new models. These scripts can be downloaded from 

http : / /projects . hepf orge .org/ sarah/ Toolbox . html 

After downloading the current version of the package, untar it using 

> tar -xzf toolbox-$Version . tar . gz 

To download and install the different packages, create and change to a working directory 
(which will later contain the installed packages) and call the configure script from there, 
e.g. 

> mkdir build 

> cd build 

> . . /configure 

The configure script will now proceed to check for the requirements of the different 
packages and download any missing files. All downloaded archives will be placed in the 
tarballs subdirectory of the directory containing the configure script from where they 
will be reused in all subsequent runs. Command line options can be used to disable specific 
packages and to point the script to custom locations of compilers and of the Mathemat ica 
kernel; a full list of those can be obtained by calling configure with the — help option. 

After configure finishes successfully, make can be called to build all configured 
packages 

> make 

After make has finished, the different packages are ready to use and can be found in 
subdirectories of the build directory, together with a suitable setup file for SSP. 

Also created by configure is a script which automates the implementation of a new 
models into the different packages 

> ./butler [model] 

where [model] is the models' SARAH name. For instance, 

> ./butler NMSSM 
resp. 
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> . /ImplementModel . sh OmegaShort/Regime-2 

can be used to install the NMSSM / resp. the left-right symmetric model discussed in 
sec. A. 2. A list of different command line options can be obtained with 

> ./butler --help 

After butler has completed, the new model is implemented into all selected packages and 
ready for use. 

8.2 Preparing a model by hand 

If more flexibility is required than is offered by butler, the model files must be generated 
and integrated into the various packages by hand. In order to demonstrate how the cor- 
responding workflow looks like, we show how the NMSSM version included in SARAH is 
implemented into the other programs 2 . 

In order to implement the MSSM and prepare for an automatic parameter scan using 
SSP, the following steps are necessary 

1. Prepare SARAH and generate all output 

(a) Define the model in SARAH and evaluate it 

in [1] : « SARAH. m; 

m 12] : « Start [ "NMSSM" ] ; 

(b) Create the source code for SPheno 
m 13] : MakeSPheno [ ] 

(c) Create model files for CalcHep and micrOMEGAs 
m [ 6 ] : MakeCHep [ ] ; 

(d) Create model files for WH IZARD 
m [51 : MakeWH I Z ARD [ ] ; 

The last three steps can also be done at once by using MakeAll [] . This 
produces the SARAH output for SPheno, CalcHep, WHIZARD, FeynArts and 
ETeX. 

2. Prepare SPheno 

(a) Create a subdirectory in your installation of SPheno3 . 1 or later and copy the 
gnerated code there 

> cd [SPheno Directory]/ 

> mkdir NMSSM 

> cp [SARAH Directory] /Output/NMSSM/EWSB/SPheno/* \ 

[SPheno Directory] /NMSSM/ 

2 For the details on the installation of the different packages we refer to the respective documentation of 
the programs. 
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(b) Compile SPheno 

> make Model=NMSSM 

3. Prepare CalcHep 

(a) Create a new project in CalcHep 

> cd [CalcHep Directory]/ 

> . /mkUsrDir NMSSM 

(b) Copy the model files from SARAH to CalcHep 

> cp [SARAH Directory] /Output /NMSSM/EWSB/CHep/ * \ 

[CalcHep Directory] /NMSSM/models/ 

(c) Start CalcHep and create an executable n.calchep program for a specific 
process using make n_calchep 

> cd NMSSM/ 

> . / calchep 

(d) Copy n_calchep and the other created files including the subdirectory aux to 
a new directory 

> cp -R results/* [Process Directory]/ 

(e) Repeat the last two steps for other processes if necessary 

4. Prepare micrOMEGAs 

(a) Create a new project in micrOMEGAs 

> cd [MicrOmegas Directory] 

> ./newProject NMSSM 

(b) Copy the CalcHep model files to the project directory 

> cp [SARAH Directory] /Output /NMSSM/EWSB/CHep/ * \ 

[MicrOmegas Directory] /NMSSM/work/models/ 

(c) Move the micrOMEGAs main file to the main directory and compile it 

> cd [MicrOmegas Directory] /NMSSM/ 

> mv work/models/CalcOmega . cpp . 

> make main=CalcOmega . cpp 

5. Prepare WHIZARD 

(a) Configure the model source 

> cd [SARAH Directory] /Output/NMSSM/EWSB/WHIZARD_Omega/ 

> ./configure WO_CONF I G= [WHIZARD Directory] /bin/ 
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(b) Compile and install the new model (with this configuration, the files will end 
up in $HOME/ . whizard) 

> make install 
8.3 Running the chain 

Now, everything is in its place, and the model can be used with the different programs: 

• Go the the SPheno directory, edit the LesHouches file and run SPheno 

> cd [SPheno Directory] /NMSSM/ 

> edit LesHouches . in . NMSSM 

> . / . . /bin/SPhenoNMSSM 

Make sure that the flags 

Block SPhenoInput # 

75 1 # HiggsBounds files 

76 1 # WHIZARD file 



are set in LesHouches . in . nMSSM in order to create the parameter file for WHIZARD 
as well as the input files for HiggsBounds. 

• Copy the spectrum file to the micrOMEGAs directory and the file for WHIZARD to 
your WHIZARD working directory 

> cp SPheno . spc .NMSSM [MicrOmegas Directory] /NMSSM/ 

> cp WHIZARD. par. NMSSM [WHIZARD Working Directory]/ 

• Run micrOMEGAs 

> cd [MicrOmegas Directory] /NMSSM/ 

> . /CalcOmega 

• Create a WHIZARD input file (e.g Input . sin) with your process setup and include 
the parameter file of SPheno by using Include ( "WHIZARD .par .NMSSM" ) . 

• Run WHIZARD 

> whizard Input. sin 

• Copy the SPheno spectrum files to the directory of your CalcHep executables and 
run n.calchep 

> cp [SPheno Directory] /bin/SPheno . spc .NMSSM [CH Process Directory]/ 

> cd [CH Process Directory] / 

> . / n_calchep 

To run n.calchep using a script you can also use the blind mode, see sec. 5.1. 



- 40 - 



9. Conclusions 



We have presented a toolchain for studying extensions of the MSSM. The model under 
consideration is implemented in SARAH which creates then source code for SPheno, thus 
facilitating the implementation of new models into SPheno in a modular way. The new 
SPheno modules calculate the mass spectrum using two-loop RGE running and one-loop 
corrections to the masses. The precision of the calculation of the mass spectrum for a large 
range of SUSY models is comparable to the known precision of the MSSM provided by the 
established spectrum generators (excluding two loop effects in the Higgs sector). Therefore, 
SARAH can be regarded as a 'spectrum-generator-generator', extending the availability of 
spectrum generators for different models dramatically. 

Apart from creating output for SPheno, SARAH is capable of creating model files 
suitable for use with CalcHep/ micrOMEGAs, FeynArts, HiggsBounds and WHIZARD. 
The masses and parameters calculated by SPheno can be used directly as input for those 
programs which then facilitate further studies of the models' phenomenology. In particular, 
the dark matter relic density predicted by the model can be simulated using CalcHep/ 
micrOMEGAs, while WHIZARD is well suited for examining collider observables involving 
the long decay chains typical for supersymmetric models. The scan of observables over 
parameter space can be automatized using the SSP package. 

As SARAH works in a very generic way, a wide range of possible models is covered. Mod- 
els with new fields and/or an extended gauge sector at the SUSY scale are supported as well 
as models with thresholds and/or gauge symmetry breaking at higher energy scales. This 
flexibility, together with the integration of different existing tools together with SARAH and 
SSP into the largely automated toolchain presented in this paper should greatly simplify 
the study of new supersymmetric models. 
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A. Input files for SARAH to generate a SPheno version with thresholds 

A.l Seesaw I 
Seesawl.m 

In the case of seesaw I, the particle content is extended by three generations of a gauge 
singlet Dr. 

Fields [ [8] ] = {vR, 3 , v, , 1 , 1}; 
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The new field is assigned a Yukawa like interaction and a mass term. In addition, there is 
an effective operator which is generated by integrating out the right handed neutrino. 

SuperPotential = 

{... ,{{1 ,Yv},{v,l ,Hu}},{{l/2,Mv},{v,v}},{{l/2,WOp},{l ,Hu,l ,Hu}}}; 

Since the field is heavy, it should not be included in the calculation of the vertices and 
masses at the SUSY scale. Therefore, we "delete" it. Note that, although deleted particles 
do not appear in the low energy spectrum and Feynman rules, they do contribute to the 
RGEs above the threshold scale. 



Delet eP art icles={v } ; 







SPheno.m 

We choose a unification of the soft-breaking mass of the scalar singlet with the other 
soft-breaking masses at the GUT scale. In addition, we want to define the values of the 
superpotential parameters at the GUT scale as input values in the LesHouches file. Fur- 
thermore, we also impose a mSugra like condition for the trilinear soft-breaking coupling. 
As the bilinear soft-breaking term does not influence the RGE running of the other param- 
eters, we can safely set it to zero. Altogether, the additional boundary conditions at the 
GUT scale are 



Bound aryHighScale={ 


{mv2, 


DIAGONAL m(T2} , 


{Mv, 


LHInput [Mv] } , 


{Yv, 


LHInput [Yv]} 


{B[Mv] , 


0}, 


{T[Yv] , 

}; 


Azero*LHInput [Yv] } 



We want to include three threshold scales: each generation of the gauge singlet should be 
integrated out at energies similar to their mass. Therefore, a good choice is 



Thresholds={ 






{Abs [MvIN[l , 1] 


,{v 


1]}}, 


{Abs[MvIN[2 ,2] 


,{v 


2]}} , 


{Abs[MvIN[3 ,3] 


,{v 


3]}} 


}; 







When thresholds are included, the boundary conditions for running up and down the 
RGEs can be defined separately for each threshold scale. To this end, it is first necessary 
to initialize the corresponding arrays 



Boundary Condi t ions Up=Tablc [ { } , { Length [Thresholds]}]; 
Boundary Condi tionsDown=T able [{} , {Length [Thresholds ] }] ; 

When a gauge singlet is integrated out, the Wilson coefficient of the effective dim-5 operator 
receives a contribution of the form 

K = -YjM N 1 Y u . (A.l) 
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When heavy superfields are integrated out, the mass splitting between the fermionic and 
scalar component is neglected. The masses are calculated individually at each threshold 
scale and saved in arrays with the names 



MassOf <> Name of Superfield 







Therefore, the contributions to the effective operator at the different scales are given by 

Boundary Conditions Down [ [ 1 ] ] = { 

{WOp[indexl , index2 ] ,WOp[indexl , indcx2 ] - Yv[l,indexl] Yv [ 1 , indcx2 ] / MassOfv [ 1 ] } 

}; 

Boundary Conditions Down [ [ 2 ] ] = { 

{WOp[indexl , index2 ] ,WOp[indexl , index2 ] - Yv[2,indcxl] Yv [ 2 , index2 ] / MassOfv [ 2 ] } 

}; 

Boundary Condi tionsDown [ [ 3 ] ] = { 

{WOp[indexl , index2 ] , - Yv [3 , index 1 ] Yv [ 3 , index2 ] / MassOfv [ 3 ] } 

}; 



A. 2 Implementation of a model with a gauge symmetry breaking scale in SARAH 

and SPheno 

This section is devoted to a detailed discussion of the implementation of a left-right super- 
symmetric model (based on that presented in [72]) into SPheno and SARAH. Since we are 
mainly interested in the implementation in SARAH, it is sufficient to simplify the model a 
bit by just choosing one threshold scale and not two as discussed in [72] . Adding the sec- 
ond threshold scale is straightforward, but would lead to some redundancy in the following 
discussion. 

A. 2.1 Summary of the model 

We here give only a short summary about the model and refer to [72] and references 
therein for more details. As adread discussed, we introduce only one threshold scale at 
which xSU(2)r x U(1)b-l gets broken 

SU(2) L x SU(2) R x U(1) B -l -> SU(2) L x U{l) Y (A.2) 

From GUT scale to SU(2) R x U{l) B -L breaking scale The MSSM particle content 
above the threshold is extended by the presence of four fields which are triplets under 
SU(2)l or SU(2)ji and which carry a B — L charge. In addition, there are two triplets 
which are uncharged under B — L. Furthermore, the right handed neutrino are part of the 
spectrum and the Higgs fields are arranged in so-called bi-doublets <3?. To get a non-trivial 
CKM matrix, we need at least two generations of $ fields. The particle content can be 
summarized as follows 
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Superfield 


generations 


SU(S) C 


SU{2) L 


SU(2) R 


U{\) B -L 


Q 


3 


3 


2 


1 


\ 

3 


Q c 


3 


3 


1 


2 


~~ 3 


L 


3 


1 


2 


1 


-1 


L c 


3 


1 


1 


2 


1 


$ 


2 


1 


2 


2 





A 

A 


1 


1 


3 


1 


2 


A 


1 


1 


3 


1 


-2 


A c 


1 


1 


1 


3 


-2 


A c 


1 


1 


1 


3 


2 


n 


1 


1 


3 


1 





n c 


1 


1 


1 


3 






The superpotential for the model reads 

W = Y Q Q<5>Q C + Y L L§L C - + fLAL + f*L c A c L c 
+ aASIA + a*A c ft c A c + aft$$ + a*fl c $<5> 

+ Ma A A + M^A C A C + M Q nn + M^fOT . (A.3) 

Below SU(2)r x U(1)b-l breaking scale Here, we are left with the MSSM plus the 
effective Weinberg operator which causes neutrino masses after EWSB. 

Boundary conditions In order to link both scales, we need the following set of boundary 
conditions: 

Y d = Y^ cos 0i - sin , Y u = —Yq cos 9 2 + Yq sin 9 2 , (A.4) 
Y e = Yl cos 6i - Yl sin 9i , y„ = -y L x cos 9 2 + Yl sin 9 2 , (A.5) 

where -R = sin(#i — 6 2 ). For the soft-trilinear couplings, Y is replaced with T in the 
expressions. For the sfermionic soft masses, we have 

2 2 2 2 / a r>\ 

m = m u c = m d c = m Qc , (A.6) 



m 2 u c 


= m % 


= m Q c 








m 2 


= m 2 e c 


= m\c 


M L 


= M R 


= M 2 . 



(A.7) 
(A.8) 

while in the Higgs sector we need the relations 

m 2 Hd = cos 2 #i(m|)ii + sin 2 #i(m|) 22 - sin 6>i cos 6>i [(m|)i 2 + (m|) 2 i] , (A.9) 
m 2 Hu = cos 2 6» 2 (m|)n + sin 2 9 2 {m%) 22 - sin 9 2 cos 9 2 [(m%) 12 + (m|) 2 i] , (A.10) 

In the gauge sector, we have to express the hypercharge coupling and the corresponding 
gaugino through 

9i = , (A.11) 

'2g 2 + 3g 2 BL 



_ 2g 2 M BL + 3g 2 BL M R 
Ml ~ 2g 2 + 3g 2 BL " < A " 12 > 
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A.2.2 Model files for SARAH 



We now discuss the input files for SARAH required for defining the model at the different 
scales. For shortness, we concentrate on the parts necessary for SPheno output and skip 
the gauge fixing terms and definition of Dirac spinors above the threshold scale. 

From GUT scale to SU(2) R x U(l) B -L breaking scale 

The vector and chiral superfields at the highest scale defining the gauge sector and particle 
content are defined by 



Gauge[[l]] = {B, U[l], bminl , gBL, False}; 

Gaugc[[2]] = {WL, SU[2], left, g2 , True } ; 

Gaugc[[3]] = {WR SU[2], right, g2 , True } ; 

Gauge[[4]] = {G, SU[3], color, g3, False}; 



Fields [[1]] = 


{{uL, dL}, 


3, 


qL , 


1/6, 


2 , 


1 , 


3}; 


Fields [ [2] ] = 


{{conj[dR], - conj[uR]}, 


3, 


qR, 


-1/6, 


1 , 


2, 


-3}; 


Fields [ [3] ] = 


{{vL, cL}, 


3, 


1L , 


-1/2, 


2 , 


1 , 


i}; 


Fields [ [4] ] = 


{{conj [cR] , — conj [vR] } , 


3, 


1R , 


1/2, 


1 , 


2, 


i}; 


Fields [ [5] ] = 


{{{HdO, Hup},{Hdm, HuO}}, 


2, 


Phi , 


0, 


2,- 


-2, 


i}; 


Fields [ [6] ] = 


{{{dcltaLp/Sqrt [2] , dcltaLpp}, 
















{dcltaLO, - dcltaLp/Sqrt [2]}} , 


1, 


deltaL , 


1, 


3 , 


1, 


i}; 


Fields [ [7] ] = 


{{{dcltaLbarm/Sqrt [2] , dcltaLbarO}, 
















{ dcltaLbarmm , — dcltaLbarm/ Sqrt [2] } } , 


1 , 


deltaLbar 


, -1, 


3 , 


1 , 


i}; 


Fields [ [8] ] = 


{{{dcltaRm/Sqrt [2] , dcltaRO}, 
















{deltaRmm, - dcltaRm/ Sqrt [2] } } , 


1 , 


deltaR , 


-1, 


1 , 


3 , 


i}; 


Fields [ [9] ] = 


{{{dcltaRbarp/Sqrt [2] , deltaRbarpp } , 
















{dcltaRbarO , - dcltaRbarp/ Sqrt [2] } } , 


1, 


deltaRbar 


, 1, 


1, 


3, 


i}; 


Fields [[10]] = 


= {{{omegaLO/Sqrt [2] , omcgaLp}, 
















{omegaLm, — omcgaLO/ Sqrt [2 ] } } , 


1, 


omegaL , 


o, 


3, 


1, 


i}; 


Fields [[11]] = 


= {{{omcgaRO/Sqrt [2] , omcgaRp}, 
















{omcgaRm, — omcgaRO/ Sqrt [2] } } , 


1, 


omegaR , 


o, 


1, 


3, 


i}; 



The superpotential reads 


SupcrPotential = { {{1, YQ} , 


{qL,qR,Phi}} , 


{{1, YL}, 


{1L ,1R , Phi}} , 


{{1, f}, 


{1L , deltaL , 1L}} , 


{{l,conj [f]} , 


{1R, deltaR ,1R}}, 


{{1 ,Mdclta} , 


{deltaL , deltaLbar }} , 


{{1 , conj [Mdclta] } , 


{deltaR , deltaRbar}} , 


{{-l/2,Mu3} , 


{Phi , Phi}} , 


{ { 1 , Momcga } , 


{ omegaL , omegaL } } , 


{ { 1 , c o n j [ Momcga ] } , 


{omegaR , omegaR} } , 


{{l,a}, 


{ deltaL , omegaL , deltaLbar } } , 


{{l,conj [a]} , 


{ deltaR , omegaR , deltaRbar}} , 


{{1 ,AlphaOm} , 


{omegaL , Phi , Phi } } , 


{{1 ,conj [AlphaOm]} , 


{omegaR, Phi , Phi}} }; 



The gauge bosons and gauginos of the right sector decompose into 



DEFINITION [RSB] [ GaugeSector] = 

{{VWR{1 ,{VWRn, 1/Sqrt[2]}, { conj [VWRm] , 1 / Sqrt [2] } } , 

{2,{VWRn, -I/Sqrt[2]}, {conj [VWRm] , I /Sqrt [2]}} , 

{3,{VWR0, 1}}}, 

{fWR,{l ,{fWRm, 1/Sqrt[2]}, {fWRp, 1 / Sqrt [2] } } , 

{2,{fWRm, -I/Sqrt[2]}, {fWRp, I / Sqrt [ 2 ] } } , 

{3,{fWR0, 1}}}}; 
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after the Q, and A fields have received their VEV 



DEFINITION [RBLSB ] [ VEVs] = 

{ {SomegaRO, {vR,l/ Sqrt [2] } , {sigmaOmR, I / Sqrt [ 2] } , {phiOmR, 1 / Sqrt [2] } } , 

{SdeltaRO, {vBL, 1 / Sqrt [2] } , {sigmaR , 1/ Sqrt [2] } , {phiR , 1 / Sqrt [2] } } , 

{SdcltaRbarO , {vBL,l/ Sqrt [2] } , {sigmaRbar , I /Sqrt [2] } ,{ phiRbar , 1 / Sqrt [ 2 ] } } } ; 



Finally, we need the rotations in the matter sector to the new mass eigenstates 



DEFINITION [RSB] [ MattcrScctor] = 




{ {{SdeltaRm, conj [ SdeltaRbarp ] } , 


{HpmlRl.ZCl}} , 


{ { SomcgaRm , conj [SomegaRp]} , 


{Hpm2Rl,ZC2}} , 


{{fB ,fWR0, FdeltaRO , FdcltaRbarO , FomegaRO} , 


{LO, ZN}}, 


{{{fWRm, FomegaRm} , {IWRp, FomegaRp}}, 


{{Lm,UM}, {Lp,UP}}}, 


{{phiR , phiRbar , phiOmR} , 


{hhR2, ZH}}, 


{{sigmaR, sigmaRbar, sigmaOmR} , 


{AhR2, ZP}}, 


{ { FvL , conj [FvR] } , 


{NO, Znu}}, 


{{SHdO.conj [SHuO]} , 


{SHOrl ,UH0}} , 


{{SHdm, conj [SHup] } , 


{SHCrl,UHC}}, 


{ {SomegaLm , conj [SomegaLp]} , 


{SOlrl ,U01}} , 


{{ SdeltaLp , conj [ SdcltaLbarm ] } , 


{SDLpRl,UDLp}} , 


{{SdcltaLpp , conj [ SdcltaLbarmm ] } , 


{SDLppRl.UDLpp}} , 


{{SdcltaLO ,conj [ SdcltaLbarO ] } , 


{SDLOrl ,UDL0}} , 


{{SdcltaRmm , conj [ SdeltaRbarpp ] } , 


{ SDRmmRl , UDRmm} } , 


{{SdeltaRO , conj [ SdcltaRbarO] } , 

}; 


{SDR0rl,UDR0}} 



Below SU(2) R x U(1)b-l breaking scale 

Below the breaking scale, only the particles and gauge groups of the MSSM survive as 
dynamical degrees of freedom. Therefore, the model file is almost identical that for the 
MSSM. We only point out the differences: 

• The superpotential contains the Weinberg operator 



SuperPotential = {...,{{ 1 ,WOp} ,{ 1 , Hu , 1 , Hu} } }; 



• The neutrinos are massive and mix among each other 

DEFINITION [EWSB] [ MatterSector]= { ...,{{FvL}, {FV, UV}}}; 



• These states form Majorana spinors 

dirac[[4]] ={Fv, FV, conj [FV] } ; 



A.2.3 Model files for SPheno output 

Two SPheno .m files, one for each scale range, are neccessary in order for SARAH to create 
the Fortran source code. While the first one is rather short, the second one includes all 
necessary boundary conditions. 
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Above SU(2)r x U(1)b-l breaking scale 

First, the regime must be nagged as an intermediate scale, and SARAH must be told its 
position in the scale hierarchy (counted from GUT to low scale). 



RegimeNr = 1 ; 
IntermediateScale 



True : 



Afterwards, we give a list with all particles which are integrated out at the threshold scale 
after gauge symmetry breaking. 



HeavyFields = {HpmlRl, ChiRl , Chalrl , hhRl, AhRl, 
FvRl, SVRrl, SH0rl[3], SHCrl [3] , 
SOlrl, SDLpRl, SDLppRl , 
SDLOrl , SDRmmRl, DR3rl , 
DLlrl , DL2rl , DL3rl , HOrl , HCrl}; 



The numbers in square brackets indicate that only the third generation and above is inte- 
grated out. 

In order to calculate the finite shifts of the gauge couplings and gaugino masses, it is 
necessary to define the gauge sector of the next scale NextGauge as well as the quantum 
number of the fields which are integrated out with respect to those gauge groups. 



NextGauge 
NextQN 





{U[l] , su 


2] , 


SU[3]} 


{ {HpmlRl, 


-1, 


1 , 


1}, 


{ChiRl , 


o, 


1 , 


1}, 


{Chalrl , 


-1, 


1 , 


1}, 


{hhRl , 


o , 


1 , 


1}, 


{AhRl, 


o, 


1 , 


1}, 


{FvRl , 


o, 


1 , 


1}, 


{SVRrl, 


o, 


1 , 


1}, 


{SHOrl , 


-1/2, 


1 , 


1}, 


{SHCrl , 


1/2, 


2 , 


1}, 


{SOlrl , 


o , 


1 , 


1}, 


{ SDLpRl, 


1 , 


1 , 


1}, 


{SDLppRl, 


2 , 


1 , 


1}, 


{SDLOrl , 


1, 


3 , 


1}, 


{SDRmmRl, 


-2, 


1 , 


1}, 


{DR3rl , 


-2, 


1 , 


1}, 


{DLlrl , 


1, 


1 , 


1}, 


{DL2rl , 


1, 


1 , 


1}, 


{DL3rl , 


1 , 


3 , 


1}, 


{HOrl , 


-1/2, 


1 , 


1}, 


{HCrl, 


1/2, 


2 , 


1} 



}; 




Finally, SARAH needs information on the vacuum conditions. There are two different 
VEVs, and therefore, we need to choose two parameters which are fixed by the tadpole 
equations. As there is no closed analytical solution for them, we give an approximation 
obtained by neglecting the soft-breaking terms. 



ParametersToSolveTadpolcs = {Mdelta, Momcga}; 
UseGivenTadpoleSolution = True; 

SubSolutionsTadpolesTree={ 

Mdelta -> - SignumMdelta acl vR/Sqrt[2] , 

Momcga -> - SignumMomega acl vBL~2/(2 Sqrt [2] vR) 

}; 
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SubSolutionsTadpolesLoop={}; 



Below SU(2)r x U(1)b-l breaking scale 



The second scale is not an intermediate scale, and hence 



RegimeNr = 2; 

IntermediateScale = False ; 







We make the following choice of free parameters of that model: to the set of standard 
mSugra parameters (mo, -M]/ 2 , Aq, tan/3, sign//), we add Bq, the superpotential parameter 
a, the signs of Mq and Ma, the two VEVs vr and vbl, and the threshold scale. These are 
defined in the blocks MINPAR and EXTPAR. 



ME\PAR= { 




{1, 


m0} , 


{2, 


ml2} , 


CO 


TanBeta} , 


{4, 


SignumMu} , 


{5, 


Azero } , 


{6, 


Bzero } , 


{7, 


SignumMomega } , 


{8, 


SignumMdelta } , 


{9, 


alnput } } ; 


EXTPAR = { 




{100, vRinput}, 


{101, vBLinput}, 


{200, TScalc}}; 



As in the, MSSM we fix the numerical values of [i and by the solving the two tadpole 
equations, 



ParametersToSolveTadpoles = {\ [Mu] ,B [\ [Mu] ] } ; 







Furthermore, we use also the common definitions for the SUSY scale, see also sec. 3.5. 



RenormalizationScaleFirstGuess = m0"2 + 4 ml2~2; 
RenormalizationScale = MSu [ 1 ] * MSu [ 6 ] ; 







We use as condition for the GUT scale 



ConditionGUTscale = {gBL=g2 , gl = g2 } ; 

The second entry is necessary because it might be that the masses of the fields breaking 
the left-right symmetry are above the unification scale which has to be handled seperatly. 
At the GUT scale we use boundary conditions motivated by minimal supergravity: all 
scalar soft-breaking masses are proportional to uiq, the gaugino masses are proportional 
to Mi/ 2, the trilinear soft-breaking couplings are given by the corresponding superpoten- 
tial parameter times Aq, and the bilinear soft-breaking couplings are set to -Bo times the 
superpotential parameter. The values for the coupling matrices /, a as well as for fj, are 
read from the LesHouches input file. 

Boundary HighScale={ 

{T[YQ] , Azcro*YQ}, 
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{T [YL] , 
{f , 

{T[f] , 
{AlphaOm, 
{T[ AlphaOm] , 
{T[a], 
{B[Mdelta] , 
{B[Momega] , 
{B[Mu3] , 
{mqL2 , 
{mqR2, 
{mlL2 , 
{mlR2, 
{mPhi2 , 
{mdeltaL2 , 
{ mdeltaLbar2 , 
{mdeltaR2 , 
{ mdeltaRbar2 , 
{momegaL2 , 
{ momegaR2 , 
{ MassB , 
{MassWL, 
{MassG , 



}; 



Azero*YL} , 

LHInput [ f ] } , 

Azero*LHInput [fm] } , 

LHInput [AlphaOm] } , 

Azero*LHInput [AlphaOm] } 

Azero*aInput} , 

Bzero* Mdclta } , 

Bzero*Momcga} , 

Bzoro*LHInput [Mu3] } , 

DIAGONAL m(T2} , 

DIAGONAL m0"2} , 

DIAGONAL mCT 2} , 

DIAGONAL m(T 2} , 

DIAGONAL m0~2}, 

m(T 2} , 

mO " 2 } , 

m(T 2} , 

mO " 2 } , 

mO " 2 } , 

mO * 2 } , 

ml2}, 

ml2}, 

m!2} 



To glue the both regimes, we need to define the appropriate boundary conditions. First, 
we initialize the arrays 



ThresholdScales = {TSCALE}; 



Boundary ConditionsUp = Table [{},{ Length [ ThresholdScales ]}] ; 
BoundaryConditionsDown = Tabic [{},{ Length [ ThresholdScales ]}] ; 



and then encode the equations eqs. (A.4)-(A.12). In order to keep the code short, we define 



STl 


= Sin 


Thetal ] ; 


CT1 


= Cos 


Thctal j ; 


ST2 


= Sin 


Thcta2 j ; 


CT2 


= Cos 


Theta2 ] ; 


ST21 


= Sin 


Theta2-Thetal ] ; 


CT21 


= Cos 


Theta2-Thetal j ; 





Using these abbreviations, the boundary conditions can be written as 



Boundary ConditionsUp [[ 1 ] ] = { 

{YQ[indexl , index2 ,1] , (Yu [ indexl , index2 ] STl + Yd [ indcxl , index2 ] ST2)/ST21 }, 

{YQ[indcxl , index2 ,2] , (Yu [ indexl , index2 ] CT1 + Yd [ indcxl , index2 ] CT2) / ST21 }, 

{YL[ indexl , index2 ,1] , (Yv [ indexl , index2 ] STl + Ye [ indcxl , indcx2 ] ST2)/ST21 }, 

{YL[indexl , indcx2 ,2] , (Yv [ indexl , index2 ] CT1 + Ye [ indcxl , index2 ] CT2)/ST21 }, 

{gBL, Sqrt[2] gl g2 /Sqrt[5 g2~2 -3 gl'2]}, 

{Yv, LHInput [Yv]} 

}; 



BoundaryConditionsDown [[ 1 ] ] = { 



{vR , vRinput } , 

{vBL, vBLinput}, 

{a , alnput } , 

{Thetal , ArcTan [ RcalPart [(( vR* AlphaOm [ 1 ,2] )/2 + Mu3 [ 1 ,2] )/Mu3[2 ,2]] ] } , 

{Theta2 , ArcTan [ RcalPart [(- (vR* AlphaOm [ 1 ,2] )/2 + Mu3 [ 1 ,2] ) /Mu3[2 ,2]]] } , 

{gl, Sqrt[5] g2 gBL/Sqrt[2 g2"2 + 3 gBL " 2] } , 

{MassB, (2 g2"2 MassB + 3 gBL '2 MassWL)/(2 g2"2 + 3 gBL ~ 2 ) } , 
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{MassWB, MassWL} , 

{Yd[indexl ,index2] , YQ[ indexl , indcx2 , 1 ] CT1 - YQ[ indexl , index2 , 2 ] ST1}, 
{Yu[indcxl ,index2] , - YQ [ indexl , index2 , 1 ] CT2 + YQ [ indexl , index2 , 2 ] ST2 } , 
{Ye[indexl ,index2] , YL [ indexl , indcx2 , 1 ] CT1 - YL [ indexl , index2 , 2 ] ST1}, 
{Yv[indexl ,index2] , - YL [ indexl , index2 , 1 ] CT2 + YL [ indexl , indcx2 , 2 ] ST2 } , 
{T[Yd] [ indexl , index2 ] , T[YQ] [ indexl , indcx2 ,1] CT1 — T[YQ] [ indexl , index2 ,2] ST1} , 
{T[Yu] [ indexl , index2 ] , -T [YQ] [ indexl , index2 , 1 ] CT2+T [YQ] [ indexl , index2 , 2 ] ST2}, 
{T[Ye] [ indexl , index2 ] , T [YL] [ indexl , index2 , 1 ] CT1 -T [YL ][ indexl , indcx2 , 2 ] ST1}, 
{mu2, mqR2}, 
{md2, mqR2}, 
{mq2, mqR2}, 
{me2, mlR2}, 
{ml2, mlR2}, 

{mHd2, CT1"2 mPhi2[l,l] + ST1"2 mPhi2[2,2] - ST1 CT1 ( mPhi2 [ 1 , 2 ] + mPhi2 [ 2 , 1 ] ) } , 
{mHu2, CT2"2 mPhi2[l,l] + ST2"2 mPhi2[2,2] - ST2 CT2(mPhi2 [ 1 , 2] + mPhi2 [ 2 , 1 ] ) } , 
{WOp, MatMul2 [ MatMul2 [Yv,InverseMatrix[f] .FortranFalse] , Transpose [Yv] , 

FortranFalse]/ vR} 

}; 



Note that Y v does not appear in any of the model files and it is therefore necessary to fix 
the dimension of that matrix by hand 



AdditionalVariablesSPheno={Yv[3 ,3] } ; 







Several parameters are restricted to be real. In addition, it is helpful to choose initialization 
values for some parameters to stabilize the numerics in the first iteration 



RealParameters = {TanBeta, vRinput , vBLinput , Thetal , Theta2 , TScale } ; 

I n i t i ali z a t i o n V al ue s = { 
{Mu3IN[l,l], (Mu3IN[l ,2] ~2 - AlphaOmIN [ 1 , 2 ] * 2 vRinput " 2/4)/Mu3IN [2 , 2] } , 
{Thetal , ArcTan [ RealPart [ - (Mu3IN [1 , 2] + AlphaOmIN [ 1 ,2] vRinput /2)/Mu3IN [2 ,2] ]] } , 
{Theta2 , ArcTan [ RealPart [(Mu3IN[l ,2] - AlphaOmIN [1 ,2] vRinput /2)/Mu3IN [2 ,2] ] ] } , 
{Mdelta, alnput * SignumMdelta* vRinput /2 }, 
{Momega, SignumMomega* ( alnput " 2* vBLinput " 2 )/( 8 Mdelta)} 

}; 



B. Example for using SSP in the MSSM: (m ,M 1 / 2 ) grid scan 

In this section we give an example on how to generate scans and plots within in the MSSM. 
The example shown here is a short version of the file Example 1 . m included in SSP. The 
goal will be varying mo and M1/2 between and 1000 GeV and thus investigating how 
different masses depend on these two parameters. 

First, it is necessary to define the location of all tools 3 in a file called Def aultSettings . m. 
The content of Def aultSettings . m should look like 

DEFAULT [SPheno] = "[ Directory ]/ SPheno3 . 1 . 4/ bin/SPhenoMSSM" ; 
DEFAULT[ SPhenoInputFilc] = " LesHouches . in .MSSM" ; 
DEFAULT[ SPhenoSpectrumFile] = " SPheno . spc .MSSM" ; 

DEFAULT [MicroOmegas] =" [Directory] / micromegas_2 .4.1 /MSSM/CalcOmega" ; 
DEFAULT[MicroOmegasInputFile] = " SPheno . spc .MSSM" ; 

3 If the environment is set up using our provided script (c.f. 8.1), the necessary file is created automatically 
in the process. 
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DEFAULT [ MicroOmegasOutputFile ] = "omg.out"; 

DEFAULT [WHEARD] = "[ Directory ]/ whizard/bin/whizard " ; 
DEFAULT [WHIZARDparFilc] = "WHIZARD. par .MSSM" ; 

DEFAULT [HiggsBounds] = "[ Directory ]/ HiggsBounds / HiggsBounds LandH effC 3 1"; 



[Directory] stands for the installation directory of the framework presented here. In 
principle, for this example it would have been sufficient to set only the information for 
SPheno. 

The properties of the desired scan have to be defined in a second input file located in 
SSP/Input which we will call ExamppleJMSSM.m. The first thing to into the file is a list 
with identifiers for the different scans, providing the possibility to perform several scans in 
a row using only one input file. Here, we will restrict ourself to just one scan 



RunScans = {ScanM0M12}; 



In order to actually perform the scans, the necessary information for creating the LesHouches 
input has to be provided, which comprises all appearing block names as well as numerical 
values for the parameters. 



DEFINITION [ScanM0M12 ] [ B 1 o c ks ] = {MODSEL, SMINPUTS , MINPAR, SPhenoInput } ; 

DEFINITION [ ScanM0M12 ] [ MODSEL] = { 
{{l},{Valuc->l}} ! 
{{6} ,{Valuc->l}} 

}; 

DEFINITION [ ScanM0M12 ] [ SMINPUTS] = { 
{{2} .{Value - > 1.166390*10 "-5}}, 
{{3},{Valuc->0.1172}}, 
{ {4}, { Value ->91. 18760}}, 
{{5} ,{Valuc->4.2}}, 
{{6},{Valuc->172.9}}, 
{{7} ,{Value->1.777}} 

}; 

DEFINITION [ ScanM0M12 ] [ MINPAR] = { 

{{1}, {Min->0,Max->1000, Steps ->10, Distribution ->LINEAR} } , 
{{2}, {Min->0,Max->1000, Steps ->10, Distribution ->LINEAR} } , 
{{3}, {Valuc->10}}, 
{{4}, {Valuc->1}}, 
{{5}, {Valuc->0}} 

}; 

DEFINITION [ScanM0M12] [ SPhenoInput]={ 
{{1} ,{Valuc->-l}}, 
{{2} ,{Valuc->l}}, 
{{ll},{Valuc->0}}, 
{{12},{Valuc->12}}, 
{{21} ,{Valuc->0}} 

}; 



Observe that all parameters are assigned constant value, except mo and M 1 / 2 which are to 
be linearly varied. 

Finally, we define the plot variables. In this example, we choose 
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• All three neutral Higgs masses for different mo as functions of Mi/ 2 

• All six selectron masses for different mo as functions of Mjy 2 

• All four neutralino masses for different Miy 2 as functions of mo 

• All three neutral Higgs masses for different Mi/ 2 as functions of mo 

• Contour plot of the light Higgs mass in the (mo, M 1 / 2 )-plane 

• Contour plot of the heavy Higgs mass in the (mo, M 1 / 2 )-plane 



DEFINITION [ScanM0M12] [ Plots] = { 

{P2D , {MINPAR[1] ,{MASS[25] ,MASS[35] ,MASS[36]}} , 
Stylel ," mO_Higgs_different_ml2 .eps"}, 

{P2D , {MINPAR[1] ,{MASS[1000011] ,MASS[200001 1] ,MASS[ 1 0000 1 3] , 
MASS[2000013] ,MASS[ 1 0000 1 5] ,MASS[20000 1 5] } } , 
Style2 , "mO_Selectrons_different_ml2.eps"}, 
{P2D , {MINPAR[2] ,{MASS[25] ,MASS[35] ,MASS[36]}} , 
Style3 , " ml2_Higgs_different_m0.eps"}, 

{P2D , {MINPAR[2] , {MASS[1000022] ,MASS[1000023] , 
MASS[1000034] ,MASS[ 1 000035] } } , 
Style4 , "ml2_Neutralinos_different_m0.eps"}, 

{P3D, {MINPAR[1] ,MINPAR[2] ,MASS[25] } , Style5 ," mO_ml2_Mass25 . eps" } , 

{P3D, {MINPAR[1] ,MINPAR[2] ,MASS[35] } , Style5 ," mO_ml2_Mass35 . eps"} 

}; 



Finally, we polish the plots by adjusting Mathematica's plot options 

Style5= {Frame— >True , Axes— >False , 

FrameLabel->{Style [Subscript ["m" ,"0"] ,16] , Style [Subscript ["M" ,"1/2"] ,16]} , 
FrameTicksStylc — > Directive [ Black , 1 4] , ContourLabels— >True } ; 



For a detailed description of those options, we refer to the manual of Mathematica. The 
final plots are shown in Fig. 4. 
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200 400 600 800 200 400 600 800 

m m 

Figure 4: From top to bottom, left to right: mass of the light Higgs in the (mo, -Mi/2)-plane, mass 
of the heavy Higgs in the (m , M 1 / 2 )-plane, masses of the scalar Higgs fields as functions of M 1 / 2 
for different values of mo, masses of the neutralinos as functions of Mij 2 for different values of mo, 
masses of the heavy Higgs in the (mo, Mi/ 2 ) plane, masses of the scalar Higgs fields as functions 
of mo for different values of M\i 2 , masses of the charged sleptons as functions of mo for different 
values of Mi/ 2 . 



- 53 - 



C. Validation of SPheno output 
C.l Mass spectrum of the MSSM 

In order to demonstrate the agreement between the "stock" SPheno MSSM and the 
SPheno version generated by SARAH, the relative mass difference 

A "JSPheno — ^SPheno — SARAH i\ 

Am = (^-1) 

m SPheno 

for two different points as functions of too resp. Mi / 2 is shown in Fig. 5 and Fig. 6 

All comparisons were done using SPheno 3.1.4 . For obtaining comparable results, the 
"stock" SPheno was slightly modified in order to compensate for the systematic differences 
shown in Tab. 1. Specifically, we restricted the calculation of the mass spectrum to two- 
loop RGEs and one- loop self-energies, thus neglecting the two- loop mass corrections in the 
Higgs sector normally included by SPheno. In addition we switched off the loop induced 
decays of neutralinos and gluinos. We used a relative precision of O(10 -4 ) to calculate the 
spectrum with both SPheno versions. 

C.2 Decay widths and branching ratios in the MSSM 

As discussed in Sec. 3, SARAH generates code for SPheno to handle all two-body decays as 
well as the fermionic three-body decays. Shown in is Tab. 4 is a comparison of the Higgs 
decays in the MSSM between SPheno and SPheno-SARAH, including all branching ratios 
which are larger than 10 -4 . A comparison of the neutralino decays into two- and three- 
body final states is given in Tab. 5. Furthermore, in Fig. 7 we show the relative differences 
between the total decay widths in SPheno "stock" vs. SPheno/SARAH as functions of mo 
for a scenario based on SPS 4. 

C.3 Masses in Seesaw type I — III 

In order to demonstrate that the automatic inclusion of thresholds in SARAH works and 
that the 1-loop boundary conditions are calculated correctly, we show a comparison for the 
seesaw type I— III between SPheno "stock" and SPheno/SARAH in Tab. 6. The model and 
SPheno .m files for this scenario are shown in appendix A.l. 

C.4 Masses in the NMSSM 

We have compared the results of a SPheno version for the NMSSM created by SARAH with 
NMSSM-Tools and found good agreement. For detailed discussion including numbers, see 
[55]. 
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Figure 5: Comparison of the masses of the Higgs fields and neutralinos between SPheno "stock" 
and SPheno/SARAH. The input values are chosen as M 1 / 2 = 500 GeV, A = GeV, tan/3 = fO, 
sign/i = 1. SPheno "stock" is based on SPheno 3.1.4 with small changes as explained in the text. 




Figure 6: Comparison of the masses of all six charged sleptons and down squarks. Input values: 
M = 250 GeV, A = -300 GeV, tan/3 = 20, sign/i = 1. 
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Figure 7: The dependence of the total decay widths of the up-squarks and the gluino on mo. 
Input values: M 1 / 2 = 500 GeV, A = GeV, tan /3 = 10, sign fj, — 1. As explained in the text, loop 
induced decays g — > \°9 were disabled in SPheno "stock". 
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hi 


h 2 


A 


r [GeV] 


2.21 • 10" 


3 2.21 • 10~ 3 


7.30- 10- 1 


7.30 • IO" 1 


1.16 


1.16 


Branching Ratios 


ss 


0.0003 


0.0003 


0.0003 


0.0003 


0.0002 


0.0002 


bb 


0.7763 


0.7760 


0.6402 


0.6403 


0.4034 


0.4032 


Hp, 


0.0005 


0.0005 


0.0004 


0.0004 


0.0002 


0.0002 


tT 


0.1287 


0.1288 


0.1051 


0.1051 


0.0662 


0.0662 


cc 


0.0380 


0.0379 


< io- 4 


< io- 4 


< 10~ 4 


< 10~ 4 


ti 


- 


- 


0.0667 


0.0668 


0.1224 


0.1223 


xixi 


- 


- 


0.0400 


0.0400 


0.2010 


0.2010 


x?x? 


- 


- 


0.0234 


0.0234 


0.0230 


0.0230 


xlxl 


- 


- 


0.0643 


0.0643 


0.0962 


0.0962 


xlxl 


- 


- 


0.0145 


0.0145 


0.0758 


0.0758 


h\h\ 


- 


- 


0.0143 


0.0143 


- 


- 


hiZ 


- 


- 


- 


- 


0.0021 


0.0021 


ejei 


- 


- 


0.0062 


0.0062 


0.0049 


0.0048 


e\e % 






0.0052 


0.0052 


< 10~ 4 


< IO" 4 


e* 2 e 2 






0.0006 


0.0006 


< IO" 4 


< IO" 4 


e* 3 e 3 






0.0006 


0.0006 


< io- 4 


< IO" 4 


e 6 ei 






0.0052 


0.0052 


0.0049 


0.0048 








0.0012 


0.0012 


< IO" 4 


< IO" 4 








0.0011 


0.0011 


< 10~ 4 


< IO" 4 








0.0011 


0.0011 


< io- 4 


< IO" 4 


77 


0.0025 


0.0025 


< io- 4 


< 10~ 4 


< io- 4 


< 10~ 4 


99 


0.0310 


0.0309 


0.0007 


0.0007 


< 10~ 4 


< 10~ 4 


ZZ* 


0.0011 


0.0010 










zz 






0.0029 


0.0029 






w-{w + y 


0.0110 


0.0110 










(w-)*w + 


0.0110 


0.0110 










w-w + 






0.0060 


0.0060 







Table 4: Left column: SPheno/ SARAH, right column: SPheno "stock". Input values: m = 
100 GeV,M 1/2 = 250 GeV, A = -100 GcV, tan £=10, signal. 
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X°2 


xl 


xl 


r [GeV] 


6.57-10" 


6 6.57 -IO" 6 


1.39 


1.39 


1.60 


1.60 


Branching Ratios 


xldd 


0.1742 


0.1742 


< 10~ 4 


< io- 4 


< io- 4 


< io- 4 


Xiss 


0.1742 


0.1742 


< 10" 4 


< io- 4 


< io~ 4 


< 10~ 4 


xibb 


0.1762 


0.1763 


< 10" 4 


< io- 4 


< io- 4 


< io- 4 


X?ee 


0.0227 


0.0227 


< IO" 4 


< io- 4 


< 10~ 4 


< io- 4 




0.0227 


0.0227 


< io- 4 


< io- 4 


< io- 4 


< 10~ 4 




0.0230 


0.0231 


< IO" 4 


< 10~ 4 


< IO" 4 


< 10~ 4 




0.1293 


0.1293 


< IO" 4 


< 10~ 4 


< 10~ 4 


< IO" 4 


Xicc 


0.1289 


0.1287 


< 10~ 4 


< IO" 4 


< IO" 4 


< IO" 4 




0.1487 


0.1486 


< IO" 4 


< 10~ 4 


< 10~ 4 


< 10~ 4 


XiW+ 






0.3137 


0.3137 


0.3519 


0.3519 


xiw- 






0.3137 


0.3137 


0.3519 


0.3519 


Xlh! 






0.0269 


0.0260 


0.0904 


0.0904 








0.1228 


0.1228 


0.1467 


0.1455 


xiz 






0.1359 


0.1359 


0.0327 


0.0327 


x\z 






0.1985 


0.1985 


0.0264 


0.0264 



Table 5: Left column: SPheno/ SARAH, right column: SPheno "stock". Input values: m = 
550 GeV,M 1/2 = 200 GeV, A = -100 GeV, tan £=10, signal. 
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Seesaw I 


Seesaw II 


Seesaw III 


particle 


SPheno 


SARAH 


SPheno 


SARAH 


SPheno 


SARAH 


-^GUT 


2.36 • 10 16 


2.36 • 10 16 


2.75 • 10 16 


2.75 • 10 16 


3.61 • 10 16 


3.64 • 10 16 


X? 


97.59 


97.61 


76.07 


76.17 


47.04 


46.92 


X% 


180.26 


180.30 


139.44 


139.51 


84.42 


84.21 


xl 


345.69 


345.93 


294.96 


294.47 


209.78 


209.56 


xl 


365.60 


365.81 


315.62 


315.22 


234.36 


234.39 


x7 


179.74 


179.78 


138.59 


138.65 


81.91 


81.70 


Xo 


366.57 


366.78 


317.00 


316.60 


236.90 


236.69 


g 


621.94 


622.17 


509.89 


509.78 


366.28 


365.81 


hi 


105.73 


105.73 


104.11 


104.11 


101.41 


101.40 


h 2 


443.91 


444.07 


384.29 


385.31 


320.91 


320.74 


A 


443.64 


443.81 


384.03 


385.31 


320.64 


320.46 


H+ 


451.26 


451.41 


392.65 


393.65 


330.69 


330.56 


e\ 


263.89 


263.85 


258.62 


258.77 


253.34 


253.38 


&2 


270.25 


270.25 


264.94 


265.09 


259.33 


259.31 


63 


270.25 


270.25 


264.97 


265.11 


259.33 


259.31 


e.4 


304.59 


304.58 


288.40 


288.46 


270.48 


270.42 


65 


304.60 


304.59 


288.41 


288.47 


270.48 


270.43 


^6 


306.38 


306.22 


290.41 


290.56 


272.41 


272.36 


d\ 


551.93 


552.14 


466.24 


466.39 


357.84 


357.51 


d~2 


592.29 


592.50 


505.04 


505.11 


396.81 


396.48 


d 3 


593.97 


594.18 


506.27 


506.35 


398.12 


397.78 


d± 


593.98 


594.18 


506.27 


506.35 


398.12 


307.79 


d§ 


614.46 


614.82 


522.97 


523.06 


409.38 


409.03 


d§ 


614.46 


614.82 


522.97 


523.07 


409.38 


409.04 


Ui 


439.39 


439.62 


367.16 


367.10 


277.73 


277.47 


U2 


594.19 


594.40 


505.87 


505.93 


396.56 


396.22 


U3 


594.20 


594.41 


505.88 


505.94 


396.56 


396.22 


U4 


609.64 


609.84 


517.06 


517.16 


401.71 


401.35 


U 5 


609.64 


609.84 


517.07 


517.16 


401.72 


401.36 


u 6 


612.32 


612.56 


530.93 


530.92 


426.07 


425.71 




292.48 


292.29 


275.68 


275.74 


256.86 


256.81 




293.79 


293.78 


276.94 


277.00 


258.14 


258.09 




293.79 


293.78 


276.94 


277.00 


258.15 


258.09 



Table 6: Comparison for Seesaw type I— III. Input values: m = 250 GeV,M!/ 2 = 250 GeV, 
A a = —0 GeV, tan ,5=10, sign/z=l and a (degenerated) Seesaw scale of 5 • 10 14 GeV was used. 
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D. Validation of the WH IZARD model files 



We compared the WH IZARD model generated by SARAH from MSSM-NoFV to the MSSM 
implementation included in WHIZARD. The result is shown in Tab. 7 for a selection of 2 — > 2 
processes and shows complete agreement between the two models. 

For the comparison, the input parameters were chosen as defined by the following 
SINDARIN code: 



mh 


_ 


1.07507347E+02 


# 


light Higgs mass 


wh 


— 


3.228E-3 


# 


light Higgs width 


mHH 


— 


4.22642020E+02 


# 


heavy Higgs mass 


mHA 


— 


4.22511683E+02 


# 


axial Higgs mass 


mHpm 


— 


4.30427526E+02 


# 


charged Higgs mass 


al_h 


— 


-1.10771279E-01 


# 


Higgs mixing angle alpha 


muJi 


— 


3.94354160E+02 


# 


Higgs mu parameter 


tanb_h 


— 


10 


# 


Higgs mixing angle tan(beta) 


msul 


— 


5.64102822E+02 


# 


u— squark mass 


msdl 


— 


5.69542907E+02 


# 


d— squark mass 


mscl 


— 


5.64114921E+02 


# 


c— squark mass 


mssl 





5.69542932E+02 


# 


s— squark mass 


mstopl 





3.65489513E+02 


# 


t— squark mass 


msbl 


— 


5.05397426E+02 


# 


b— squark mass 


msu2 


_ 


5.46608780E+02 


# 


u— squark mass 


msd2 


— 


5.46345296E+02 


# 


d— squark mass 


msc2 


— 


5.46592401E+02 


# 


c— squark mass 


mss2 


= 


5.46340940E+02 


# 


s— squark mass 


mstop2 




5.84545812E+02 


# 


t— squark mass 


msb2 




5.45062834E+02 


# 


b— squark mass 


msel 


_ 


1.89600608E+02 


# 


selectronl mass 


msne 





1.72254490E+02 


# 


elect ron —sneutrino mass 


msmul 


= 


1.89622816E+02 


# 


smuonl mass 


msnmu 




1.72247165E+02 


# 


muon— sneutrino mass 


mstaul 




1.08017023E+02 


# 


staul mass 


msntau 


_ 


1.70180182E+02 


# 


tau— sneutrino mass 


mse2 


_ 


1.25413287E+02 


# 


selectron2 mass 


msmu2 


= 


1.25349593E+02 


# 


smuon2 mass 


mstau2 


= 


1.94644276E+02 


# 


stau2 mass 


mgg 




6.14828189E+02 


# 


gluino mass 


mchl 


= 


1.83712010E+02 


# 


charginol mass (signed) 


mch2 




4.13907598E+02 


# 


chargino2 mass (signed) 


mneul 




9.79820056E+01 


# 


neutralinol mass (signed) 


mneu2 




1.83937651E+02 


# 


neutralino2 mass (signed) 


mneu3 




-3.98749572E+02 


# 


neutralino3 mass (signed) 


mneu4 




4.12379874E+02 


# 


neutralino4 mass (signed) 


mt.ll 




5.64067874E-01 


# 


Rc[R.st (1,1)] 


mt_12 




8.25728426E-01 


# 


Rc[R_st (1,2)] 


mt_21 




-8.25728426E-01 


# 


Rc[R_st (2,1)] 


mt_22 




5.64067874E-01 


# 


Rc[R_st (2,2)] 


mb_ll 




9.39288883E-01 


# 


Rc[R_sb (1 ,1)] 


mb_12 




3.43127372E-01 


# 


Rc[R_sb (1 ,2)] 


mb_21 




-3.43127372E-01 


# 


Rc [ R_sb (2,1)] 


mb_22 




9.39288883E-01 


# 


Rc[R_sb (2 ,2)] 


mLll = 




3.14734480E-01 


# 


Re[R_sta (1 ,1)] 


ml_12 = 




9.49179755E-01 


# 


Re[R_sta (1 ,2)] 
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ml 


2 1 




— Q 4Q1 7Q7^^V— 01 


44 
it 


Re 




_sta (2 


ml 


22 




3 1 4 7 34 4 8 OF,— D1 


44 
TT 


Re 


R 


_sta (2 




1 1 




— Q 8fi90344fiF— 01 


44 
7T 


Re 


N ( 


1 


1)] 




12 




5 57184333F— 02 


44 
TT 


Re 


N ( 


1 


2)] 




13 




— 1 494&& c i1 *3F— 01 


44 

it 


Re 


N ( 


1 


3)] 




14 




5 1 1 27Qfi c S1 F— f)2 


44 
TT 


Re 


N ( 


1 


4)] 




21 




— Q 81 ^1 Q*37QF— 09 


44 
it 


Re 


N ( 


2 


1)] 




22 




— Q ^08^7^77^—01 


44 
it 


Re 


N ( 


2 


2 ) ] 




23 




2 C S704(S44QE— 01 


44 
TT 


Re 


N ( 


2 


3)] 


mn 


.24 


= 


-1.41997141E-01 


# 


Re 


N( 


2 


4)] 


mn 


.31 


= 


-5.87383665E-02 


# 


Re 


N( 


3 


1)] 




32 




8 K9 c »fiQ9 ^F— 09 


it 


Re 


N ( 


3 


2)] 




33 




QR119QfifiF— 01 


44 
it 


Re 


N ( 


3 


3)1 




34 




7 1 nnfi78 1 3F,— m 

1 . 1 Vl U \J 1 O 1 U i-J W _L 


44 
it 


Re 


N ( 


3 


4)] 


mn 


.41 


= 


-1.13743529E-01 


# 


Re 


N( 


4 


1)] 


mn 


.42 


= 


2.91522006E-01 


# 


Re 


N( 


4 


2)] 




43 




fi ^ n 1 9 7R0F- f)1 


44 
it 


Re 


N ( 


4 


3)] 




11 




— f\ 877RQ17QF— 01 


it 


Re 


N ( 


4 


4)] 


mu 


.11 




-9.34453143E-01 


# 

TT 


Re 


U( 


1 


1)] 


mu 


.12 




3.56086119E-01 


# 


Re 


Ui 


1 


2)] 


mu 


.21 




3.56086119E-01 


# 


Re 


U( 


2 


1)] 


mu 


.22 




9.34453143E-01 


# 


Re 


U( 


2 


2)] 


mv 


.11 




-9.81098181E-01 


# 


Re 


V( 


1 


!)] 


mv 


.12 




1.93510619E-01 


# 


Re 


V( 


1 


2)] 


mv 


.21 




1.93510619E-01 


# 


Re 


V( 


2 


1)] 


mv 


.22 




9.81098181E-01 


# 


Re 


V( 


2 


2 )] 



,1)] 

,2)] 



The intergration was performed for a CMS energy of 2 GeV, with cuts and integration pa- 
rameters chosen as exemplified by the following sample input (with the obvious replacement 
of 11/2 and F1/F2 by the actual particle identifiers): 

model = MSSM 

include (" par_MSSM . txt " ) 

process scatter = 11,12 => F1,F2 

compile 

sqrts = 2000 GeV 

cuts = all E >= 10 GeV [Fl] and 
all E >= 10 GeV [F2] and 
all 0.1<Thcta<2 [Fl] and 
all 0.1<Thcta<2 [F2] 

integrate ( scatter ) {iterations = 10:5000,10:25000,10:25000:""} 
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Process 


r n i 

cr 5 [fb] 


5 S [fb] 


a w [fbj 


Sw [fb] 


a r n i 

A [fb] 


~o ~n 

X1X1 


-)• ZZ 


1.227 x 10" 1 


9.87 


x 


10" 6 


1.227 x 10" 1 


1.35 


x 


10" 5 


3.28 x 10- 5 


~n ~ n 

X1X1 


~n — n 

-> X1X1 


2.538 x 10" 2 


1.34 


x 


10" 6 


2.538 x 10" 2 


3.65 


x 


10" 6 


-2.99 x 10" 6 


~n ~n 
X1X1 


~n ~n 
-> X2X2 


2.219 x 10" 1 


1.14 


x 


10" 5 


2.219 x 10" 1 


1.03 


x 


10" 5 


3.4 x 10" 6 


~n ~n 

X1X1 


~n - n 

^XsXa 


1.762 


1.24 


X 


10" 4 


1.762 


1.24 


x 


10" 4 


-2.7 x 10" 4 


~n ~n 

X1X1 


-> X4X4 


1.105 


8.23 


X 


10" 5 


1.104 


6.36 


x 


10" 5 


1.37 x 10" 4 


~n ~n 

X1X1 




3.913 x 10" 1 


3.46 


X 


10" 5 


3.913 x 10" 1 


2.36 


x 


10" 5 


1.43 x 10" 5 


~n ~n 

X1X1 


-)• if A 


1.843 


2.10 


X 


10" 4 


1.841 


1.68 


x 


10" 4 


1.38 x 10" 3 


~n ~n 
X1X1 


->• x^xr 


1.069 x 10" 1 


6.69 


X 


10" 6 


1.068 x 10" 1 


1.10 


x 


10" 5 


4.44 x 10" 5 


x?x? 


-> x^x 2 


1.542 


1.69 


X 


10" 4 


1.541 


1.34 


x 


10" 4 


1.19 x 10" 3 


X1X1 


^x?x° 


1.504 x 10" 1 


1.32 


X 


10" 5 


1.504 x 10" 1 


2.04 


x 


10" 5 


-6.14 x 10" 5 


~n ~n 

X1X1 


-)• W~W + 


1.57 


1.73 


X 


10" 4 


1.57 


1.55 


x 


10" 4 


-5.73 x 10" 5 


~n ~n 

X1X1 


— > (itJ 


6.836 x 10" 1 


8.44 


X 


10" 5 


6.838 x 10" 1 


5.21 


x 


10" 5 


-2.12 x 10" 4 


~n ~n 

X1X1 


-)• 


1.234 x 10 1 


1.18 


X 


10 _ci 


1.234 x 10 1 


8.07 


x 


10" 4 


6.67 x 10" 4 


~n ~n 

X1X1 


66 


7.878 x 10" 1 


9.42 


X 


10" 5 


7.879 x 10" 1 


8.46 


x 


10" 5 


-3.97 x 10" 5 


~n ~n 

X1X1 


-> tit* 


9.561 


1.02 


X 


10" 3 


9.563 


8.33 


x 


10" 4 


-2.5 x 10" 3 


~n ~n 

X1X1 


-»• ^2 


5.421 x 10" 1 


3.56 


X 


10" 5 


5.421 x 10" 1 


4.61 


x 


10" 5 


6.33 x 10" 5 


~n ~o 
X1X1 


-> ii*2 


2.096 


1.88 


X 


10" 4 


2.098 


2.22 


x 


10" 4 


-1.92 x 10" 3 


~n ~n 

X1X1 




2.097 


1.70 


X 


10" 4 


2.097 


1.63 


x 


10" 4 


-3.18 x 10" 5 


~n ~n 

X1X1 


— >■ Jljj 


1.844 x 10" 1 


1.20 


X 


10" 5 


1.843 x 10" 1 


1.58 


x 


10" 5 


6.5 x 10" 5 


~n ~n 

X1X1 


->■ J 2 4 


7.987 x 10" 1 


7.22 


X 


10" 5 


7.986 x 10" 1 


7.23 


x 


10" 5 


1.01 x 10" 4 


x?x? 


-> 616^ 


2.059 x 10" 1 


1.54 


X 


10" 5 


2.058 x 10" 1 


2.20 


x 


10" 5 


3.32 x 10" 6 


X1X1 


-»• 6 2 6J 


2.057 x 10" 1 


2.08 


X 


10" 5 


2.06 x 10" 1 


2.82 


x 


10" 5 


-2.11 x 10" 4 


~n ~n 
X1X1 


-)• l/ e ^ e 


2.209 


2.56 


X 


10 4 


2.21 


1.65 


x 


10 4 


-1.26 x 10" 3 


~n ~n 

X1X1 




2.211 


1.59 


X 


10" 4 


2.21 


1.36 


x 


10" 4 


6.83 x 10" 4 


~n ~n 

X1X1 


-)• eie^ 


1.074 


7.01 


X 


10" 5 


1.073 


7.70 


x 


10" 5 


2.87 x 10" 4 


~n ~n 

X1X1 


->■ e 2 e2 


2.605 x 10 1 


2.79 


X 


10" 3 


2.604 x 10 1 


1.59 


x 


10" 3 


9.78 x 10" 3 


~o ~n 
X1X1 


-> nrf 


2.148 x 10 1 


1.47 


X 


lO" 3 


2.148 x 10 1 


2.17 


x 


10" 3 


5.8 x 10" 3 


~n ~n 

X1X1 


-> T 2 f 2 


4.171 x 10" 1 


2.00 


X 


10" 4 


4.169 x 10" 1 


3.14 


x 


10" 5 


1.9 x 10" 4 


~n ~n 

X1X1 


-)• fif 2 


2.64 


2.34 


X 


10" 4 


2.64 


2.11 


x 


10" 4 


5.81 x 10" 4 


x?x? 


-> 7^2 rjf 


2.639 


1.95 


X 


10" 4 


2.64 


2.05 


x 


10" 4 


-7.78 x 10" 4 


x?x? 


-)• H~H+ 


1.14 


8.12 


X 


10" 5 


1.14 


1.00 


x 


10" 4 


-4.47 x 10" 4 


x?x? 


-)• /i/i 


7.951 x 10" 2 


8.62 


X 


10" 6 


7.948 x 10" 2 


7.51 


x 


10" 6 


2.98 x 10" 5 


~n ~n 

X1X1 


-)• iiii 


5.380 x 10" 2 


9.95 


X 


10" 6 


5.380 x 10" 2 


5.76 


x 


10" 6 


3.25 x 10" 6 


~n ~o 
X1X1 


/iF 


2.643 x 10" 1 


2.19 


X 


10" 5 


2.643 x 10" 1 


1.79 


x 


10" 5 


8.54 x 10" 5 


xSx? 


-> A A 


5.487 x 10" 2 


6.71 


X 


10" 6 


5.492 x 10" 2 


1.10 


x 


10" 5 


-5.63 x 10" 5 


ee — >■ 


ZZ 


2.325 


2.42 


X 


10- 4 


2.325 


1.82 


x 


10- 4 


-1.77 x 10" 4 


ee — > 


77 


1.636 x 10 1 


1.52 


X 


10" 3 


1.637 x 10 1 


2.81 


x 


10" 3 


-1.09 x 10" 2 


ee — > 


— — * 

eiej 


5.823 


4.33 


X 


10- 4 


5.825 


4.47 


x 


10- 4 


-1.97 x 10" 3 


ee — > 


e 2 e!; 


4.955 


6.10 


X 


10- 4 


4.953 


5.79 


x 


10- 4 


1.74 x 10" 3 


ee — > 


X?X? 


1.143 x 10 1 


7.15 


X 


10- 4 


1.143 x 10 1 


1.01 


X 


10- 3 


-9.57 x 10" 4 



(continued on next page) 
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Tab. 7: continued 



Process 

ee 



xgxg 



ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 
ee - 

TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 



X° 3 X°3 

xlxl 

hA° 
HA° 

xtxi 

X2X2 

X°lX°2 

w~w + 

dd 
ti 
bb 

ht\ 

t2t* 2 
t\t* 2 

d 2 d\ 

hk 
T2T2 

f\T2 
T2T\ 

lZ 

• ZZ 

■ 77 

■xlxl 

■X°2X°2 

■x° 3 xl 
■xlxl 

■ hA° 

• HA° 

■ xtxi 

■ X2X2 
■xlxl 

■ w~w+ 

• dd 



<JS [fb] 

8.503 
1.488 x lO" 3 
4.235 x 10" 2 
1.932 x 10" 4 

1.351 

1.658 x 10 1 
9.619 
4.361 

2.656 x 10 1 

7.563 
1.428 x 10 1 

7.561 

3.695 

3.091 

7.906 x 10" 1 

7.907 x 10" 1 
3.087 

6.913 x 10" 1 
3.482 x 10" 1 
3.482 x 10" 1 

2.038 

3.382 

3.818 
1.56 x 10" 1 
1.56 x 10" 1 
1.205 x 10 1 

3.148 

2.325 
1.637 x 10 1 
1.149 x 10 1 

8.577 
1.805 x 10" 2 
8.165 x 10" 2 
6.969 x 10" 4 
9.225 x 10" 1 

1.659 x 10 1 
9.41 
4.329 

2.66 x 10 1 
7.56 



5s [fb] 
6.86 x 10" 4 

1.37 x 10" 7 

2.38 x 10" 6 

1.41 x 10" 8 
1.06 x 10" 4 

1.25 x 10" 3 

7.71 x 10" 4 

5.20 x 10" 4 
5.08 x 10" 3 
8.01 x 10" 4 

1.58 x 10" 4 
7.10 x 10" 4 
1.15 x 10" 3 

5.84 x 10" 5 
4.15 x 10" 6 
2.29 x 10" 5 

1.42 x 10" 5 

3.26 x 10" 6 
7.06 x 10" 5 

5.59 x 10" 5 
2.50 x 10" 5 
1.47 x 10" 5 
2.15 x 10" 4 

2.21 x 10" 6 
2.12 x 10" 6 

1.72 x 10" 3 

1.85 x 10" 4 

2.39 x 10" 4 

1.27 x 10" 3 

1.43 x 10" 3 
4.95 x 10" 4 
1.42 x 10" 6 
6.90 x 10" 6 
7.26 x 10" 8 

5.97 x 10" 5 
1.03 x 10" 3 
5.50 x 10" 4 

3.15 x 10" 4 

1.98 x 10" 3 

7.16 x 10" 4 



<?w [fb] 

8.503 
1.488 x 10" 3 
4.234 x 10" 2 
1.932 x 10" 4 

1.351 

1.658 x 10 1 
9.617 
4.36 

2.657 x 10 1 

7.558 
1.428 x 10 1 

7.559 

3.693 

3.087 

7.906 x 10" 1 

7.907 x 10" 1 
3.086 

6.913 x 10" 1 
3.484 x 10" 1 
3.484 x 10" 1 

2.038 

3.383 

3.816 
1.56 x 10" 1 
1.56 x 10" 1 
1.204 x 10 1 

3.147 

2.325 
1.637 x 10 1 
1.148 x 10 1 

8.576 
1.805 x 10" 2 
8.161 x 10" 2 
6.967 x 10" 4 
9.223 x 10" 1 

1.659 x 10 1 
9.412 
4.326 

2.66 x 10 1 
7.56 



Sw [fb] 
6.09 x 10" 4 

1.01 x 10" 7 
2.45 x 10" 6 
3.38 x 10" 9 

2.44 x 10" 5 

1.12 x 10" 3 
5.87 x 10" 4 

2.45 x 10" 4 
1.50 x 10" 3 

4.25 x 10" 4 
2.03 x 10" 3 
4.15 x 10" 4 

1.48 x 10" 4 

3.37 x 10" 4 

2.13 x 10" 6 
6.09 x 10" 6 
1.09 x 10" 4 
2.64 x 10" 6 
1.21 x 10" 6 

9.92 x 10" 6 

1.49 x 10" 5 
1.07 x 10" 3 
1.73 x 10" 4 
3.58 x 10" 6 

1.38 x 10" 6 
8.95 x 10" 3 
2.94 x 10" 5 
1.67 x 10" 4 

1.26 x 10~ 3 
7.34 x 10" 4 
6.90 x 10" 4 
1.32 x 10" 6 
7.69 x 10" 6 

1.02 x 10" 7 
9.79 x 10" 5 
1.23 x 10" 3 
1.01 x 10" 3 
3.45 x 10" 4 

1.93 x 10" 3 
4.47 x 10" 4 



A [fb] 
-9.71 x 10- 5 
1.28 x 10" 7 
8.23 x 10" 7 
1.21 x 10" 8 
-2.41 x 10" 4 
6.14 x 10" 3 
2.26 x 10" 3 
1.14 x 10" 3 
-1.21 x 10" 2 
5.42 x 10" 3 
-1.05 x 10" 3 
2.13 x 10" 3 
1.92 x 10" 3 
3.88 x 10" 3 
-1.3 x 10" 5 
6.73 x 10" 5 
5.13 x 10" 4 
-6.95 x 10" 7 
-2.33 x 10" 4 
-2.07 x 10" 4 
-2.14 x 10" 5 
-1.55 x 10" 3 

1.62 x 10" 3 
3.54 x 10" 7 
-5.22 x 10" 6 

1.28 x 10" 2 

1.63 x 10" 3 
-4.3 x 10" 4 
5.26 x 10" 3 
4.91 x 10" 3 
1.19 x 10" 3 

-1.36 x 10" 6 
3.49 x 10" 5 
1.57 x 10" 7 
1.57 x 10" 4 

-2.02 x 10" 3 
-1.6 x 10" 3 
2.57 x 10" 3 

-4.75 x 10" 3 
5.37 x 10" 4 



(continued on next page) 
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Tab. 7: continued 



Process 
tt 



TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 

dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 
dd 



bb 

t 2 t\ 
tit* 2 

t 2 t\ 

did\ 

d 2 d 2 

bill 

b 2 b\ 

VeV e 

e 2 e* 2 
t\t{ 

T 2 f 2 
T\f 2 
f 2 f{ 

1 Z 

H-H+ 

hh 

HH 

hH 

A°A° 

ZZ 

77 

xSx? 

X° 2 X° 2 

xlxi 

hA° 
HA° 

xtxi 

X 2 X 2 

w~w + 
dd 
tt 
bb 



os [fb] 
1.428 x 10 1 

8.38 

3.804 

3.187 
9.566 x 10" 1 
9.572 x 10" 1 

3.087 
6.913 x 10" 1 
3.484 x 10" 1 
3.48 x 10" 1 

2.036 

2.038 

4.052 

3.444 

3.454 

4.131 

1.426 

1.425 
1.205 x 10 1 

2.964 
3.539 x 10" 5 
1.317 x 10" 3 
5.755 x 10" 4 
1.256 x 10" 3 

3.043 
6.740 x 10" 2 
3.762 x 10" 2 

1.547 
9.423 x 10" 4 
2.998 x 10" 2 
9.509 x 10" 5 
6.653 x 10" 1 

4.149 

1.532 
1.732 x 10" 1 

8.039 
5.63 x 10 3 
4.218 x 10 2 
4.179 x 10 2 
1.493 x 10 2 



5s [fb] 

5.20 x 10" 4 
5.27 x 10" 4 

1.07 x 10" 4 

1.26 x 10" 5 
1.22 x 10" 4 
1.80 x 10" 6 
7.22 x 10" 5 
2.90 x 10" 5 

7.64 x 10" 7 
1.70 x 10" 4 
4.42 x 10" 4 
7.67 x 10" 5 
5.98 x 10" 5 

3.51 x 10" 4 
1.92 x 10" 4 

2.21 x 10" 4 

1.08 x 10" 4 

9.58 x 10" 5 
3.18 x 10" 3 
2.41 x 10" 4 
1.86 x 10" 9 
1.32 x 10" 7 
5.17 x 10" 8 

1.65 x 10" 7 
2.57 x 10" 4 
1.01 x 10" 5 
2.54 x 10" 6 
1.24 x 10" 4 

7.59 x 10" 8 

2.27 x 10" 6 
7.01 x 10" 9 
3.82 x 10" 5 
2.53 x 10" 4 
9.89 x 10" 5 
1.20 x 10" 5 
8.13 x 10" 4 

5.52 x 10" 1 
5.85 x 10" 2 
3.76 x 10" 2 
2.27 x 10" 3 



<?w [fb] 
1.428 x 10 1 

8.377 

3.803 

3.185 
9.564 x 10" 1 
9.569 x 10" 1 

3.085 
6.912 x 10" 1 
3.485 x 10" 1 
3.484 x 10" 1 

2.038 

2.038 

4.052 

3.447 

3.453 

4.131 

1.426 

1.425 
1.206 x 10 1 

2.963 
3.540 x 10" 5 
1.316 x 10" 3 
5.757 x 10" 4 
1.257 x 10" 3 

3.042 
6.733 x 10" 2 
3.763 x 10" 2 

1.547 
9.424 x 10" 4 
3.000 x 10" 2 
9.514 x 10" 5 
6.653 x 10" 1 

4.148 

1.532 
1.731 x 10" 1 

8.037 
5.63 x 10 3 
4.222 x 10 2 
4.176 x 10 2 
1.493 x 10 2 



Sw [fb] 
1.45 x 10" 4 

4.90 x 10" 4 
5.25 x 10" 4 
7.57 x 10" 4 

8.73 x 10" 5 
4.59 x 10" 5 
3.80 x 10" 5 
4.63 x 10" 6 
1.65 x 10" 5 
1.28 x 10" 6 

8.12 x 10" 5 
3.27 x 10" 5 

2.06 x 10" 4 

1.13 x 10" 5 
1.62 x 10" 4 
1.98 x 10" 4 
1.04 x 10" 4 
8.42 x 10" 5 
2.10 x 10" 3 
3.02 x 10" 4 
2.75 x 10" 9 
1.56 x 10" 7 
6.97 x 10" 8 

3.07 x 10" 7 
1.50 x 10" 4 
9.16 x 10" 6 
3.39 x 10" 6 

1.91 x 10" 4 
6.25 x 10" 8 
3.31 x 10" 6 
1.83 x 10" 8 
1.87 x 10" 6 

2.74 x 10" 4 
1.02 x 10" 4 
1.33 x 10" 5 
8.19 x 10" 4 
3.47 x 10" 1 
5.25 x 10" 3 
8.24 x 10" 2 
6.30 x 10" 3 



A [fb] 
-3.12 x 10" 3 

2.96 x 10" 3 
4.57 x 10" 4 
2.04 x 10" 3 
1.72 x 10" 4 
2.8 x 10" 4 
1.29 x 10" 3 

6.97 x 10" 5 
-6.59 x 10" 5 
-4.05 x 10" 4 
-1.22 x 10" 3 
-2.33 x 10" 4 
-1.19 x 10" 4 
-2.49 x 10" 3 

1.39 x 10" 3 
1.63 x 10" 4 
-3.02 x 10" 5 
2.39 x 10" 4 
-4.33 x 10" 3 
3.63 x 10" 4 
-1.12 x 10" 8 
6.19 x 10" 7 
-1.80 x 10" 7 
-5.60 x 10" 7 
6.63 x 10" 4 
6.95 x 10" 5 
-9.14 x 10" 6 
5.74 x 10" 4 
-4.87 x 10" 8 
-1.77 x 10" 5 
-4.92 x 10" 8 
-3.64 x 10" 6 
8.91 x 10" 4 
2.56 x 10" 4 
3.95 x 10" 5 
1.94 x 10" 3 
-1.87 x 10" 2 
-4.27 x 10" 1 
3.12 x 10" 1 
2.93 x 10" 2 



(continued on next page) 
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Tab. 7: continued 



Process 


cr 5 [fb] 




6 S [fb] 


r n 1 

(JVK [fbj 


Sw [fb] 


A Ffl 1 

A [fb] 


dd - 


» ?2?2 


9.922 x 10 1 


1 


34 


X 


10~ 3 


9.922 x 10 1 


4.09 


x 


10" 4 


-2.94 x 10" 3 


dd- 




3.893 x 10" 1 


4 


45 


X 


10" 6 


3.893 x 10" 1 


1.49 


x 


10" 4 


-3.14 x 10" 6 


dd- 


> tati 


3.893 x 10" 1 


1 


91 


X 


10" 5 


3.892 x 10" 1 


2.07 


X 


10" 6 


2.72 x 10" 5 


dd- 


> die£j[ 


3.286 x 10 2 


1 


63 


X 


10~ 2 


3.286 x 10 2 


2.04 


X 


^ „ 9 

10~ 2 


4.91 x 10" 2 


dd- 


> d 2 d£ 


3.973 x 10 2 


2 


08 


X 


10~ 2 


3.973 x 10 2 


2.62 


X 


10~ 2 


2.6 x 10" 2 


dd- 


> bi&2 


1.713 x 10" 1 


1 


58 


X 


10" 6 


1.715 x 10" 1 


1.97 


X 


10" 6 


-2.02 x 10" 4 


dd- 


> 5 2 6J 


1.715 x 10" 1 


2 


92 


X 


10" 6 


1.715 x 10" 1 


6.87 


X 


10" 6 


-1.12 x 10" 5 


dd- 




1.003 


4 


84 


X 


10" 5 


1.002 


6.54 


X 


10" 4 


1.22 x 10" 3 


dd- 




1.003 


1 


03 


X 


10" 4 


1.003 


1.63 


X 


10" 4 


5.07 x 10" 4 


dd- 


> eiej 


5.843 x 10" 1 


4 


91 


X 


10" 6 


5.838 x 10" 1 


1.38 


X 


10" 5 


4.89 x 10" 4 


dd- 


> e 2 e?. 


1.28 x 10" 1 


4 


74 


X 


10" 6 


1.28 x 10" 1 


3.03 


X 


10" 6 


1.72 x 10" 5 


dd- 


> nf* 


9.762 x 10" 2 


3 


05 


X 


10" 5 


9.772 x 10" 2 


3.57 


X 


10" 7 


-1.01 x 10" 4 


dd- 


> T2T2 


4.616 x 10" 1 


4 


36 


X 


10" 5 


4.617 x 10" 1 


7.55 


X 


10" 5 


-1.35 x 10" 4 


dd- 


> flT2 


7.682 x 10" 2 


9 


73 


X 


10~ 7 


7.682 x 10" 2 


4.04 


X 


10" 7 


-3.85 x 10" 6 


dd- 




7.679 x 10" 2 


7.09 


X 


10" 6 


7.683 x 10" 2 


2.60 


X 


10" 6 


-4.29 x 10" 5 


dd- 


> 7 Z 


6.594 x 10" 1 


5 


08 


X 


10" 5 


6.592 x 10" 1 


4.78 


X 


10" 5 


1.83 x 10" 4 


dd- 


> H~H + 


4.539 x 10" 1 


2 


27 


X 


10" 5 


4.539 x 10" 1 


7.42 


X 


10" 6 


-6.51 x 10" 5 


dd- 


> 59 


4.938 x 10 2 


2 


83 


X 


10~ 2 


4.94 x 10^ 


4.32 


X 


10" 2 


-2.43 x 10" 1 


dd- 


> 95 


1.118 x 10 3 


1 


03 


X 


10" 1 


1.118 x 10 3 


7.46 


X 


10" 2 


-1.83 x 10" 1 


77 " 


■> x^xr 


3.402 x 10 1 


4 


19 


X 


10~ 3 


3.401 x 10 1 


2.50 


X 


10" 3 


6. x 10" 3 


77 " 


■> W~W + 


3.22 x 10 2 


2 


21 


X 


10" 2 


3.22 x 10 2 


3.20 


X 


10" 2 


5.72 x 10" 2 


77 " 


■* x 2 x 2 


3.696 x 10 1 


2 


23 


X 


10~ 3 


3.696 x 10 1 


2.88 


X 


10" 3 


4.95 x 10" 3 


77 " 


■> dd 


1.212 


1 


75 


X 


10" 4 


1.213 


2.79 


X 


10" 4 


-6.44 x 10" 4 


77 " 


■> it 


2.008 x 10 1 


3 


19 


X 


10~ 3 


2.009 x 10 1 


2.72 


X 


10" 3 


-8.75 x 10" 3 


77 " 


■> 66 


1.213 


1 


69 


X 


10" 4 


1.213 


9.21 


X 


10" 5 


1.13 x 10" 4 


77 " 


■> *x*i 


6.059 


6 


91 


X 


10" 4 


6.058 


4.91 


X 


10" 4 


7.36 x 10" 4 


77 " 


■> i 2 i 2 


3.777 


3 


17 


X 


10 


3.778 


3.37 


X 


10 


-3.11 x 10" 4 


77 " 


->■ did* 


2.44 X 10 


1 


57 


X 


10" 5 


2.44 x 10 


1.55 


X 


10" 5 


5.97 x 10 


77 " 


■> d 2 d*. 


2.57 x 10" 1 


1 


99 


X 


10" 5 


2.569 x 10" 1 


1.54 


X 


10" 5 


7.42 x 10" 5 


77 " 


■> eicj 


1.319 x 10 1 


8 


09 


X 


10- 4 


1.319 x 10 1 


9.91 


X 


10- 4 


1.09 x 10" 3 


77 " 


■> e 2 e^ 


1.391 x 10 1 


7 


59 


X 


10- 4 


1.391 x 10 1 


1.32 


X 


10" 3 


-2.98 x 10" 4 


77 " 


■> riff 


1.406 x 10 1 


1 


52 


X 


10" 3 


1.406 x 10 1 


1.13 


X 


10" 3 


-6.26 x 10" 3 


77 " 


■> T 2 T 2 


1.312 x 10 1 


6 


91 


X 


10- 4 


1.312 x 10 1 


1.17 


X 


10" 3 


1.12 x 10" 3 


77 " 




8.983 


1 


07 


X 


10" 3 


8.978 


8.81 


X 


10- 4 


4.25 x 10" 3 



Table 7: Comparision between the WHIZARD stock MSSM and the version generated 
by SARAH for several for various 2 — > 2 processes. The calculated values for the cross 
sections are shown together with the respective absolute integration errors. See the text 
for details. 
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